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SUMMARY 


-nA  linearized  treatment  of  three-dimensional  oscillatory  flow  in 
supercritical  nozzles  has-been  performed  for  the  two  cases  where  the 
steady-state  flow  is  axisymmetric  or  two-dimensional.  In  the  axi- 
syaaetric  cose,  perturbation  series  have-been  employed  to  study  the  non¬ 
linear  oscillations.  In  these  analyses,  variables  have-been  separated, 
reducing  '.'Se  system  of  partial  differential  equations  to  a  system  of 
ordinary  differential  equations.  The  variables  describing  the  transverse 
dependencies  of  the  flow  properties  are  governed  by  well-known  differ¬ 
ential  equations  (for  example,  Bessel’s  equation  is  obtained).  The 
axial  dependencies,  on  the  other  hand,  are  governed  by  differential 
equations  which  must  be  solved  by  numerical  means.  The  nozzle  admittance 
coefficients  are  related  tc  the  axial  dependencies  of  the  flow  properties. 
Certain  techniques  have-been  applied  to  reduce  the  order  of  the  differ¬ 
ential  equations  for  the  purpose  of  easier  calculation  of  the  admittance 
coefficients. 

These  admittance  coefficients  are  to  be  used  in  the  boundary  condition 
applied  at  the  exit  of  the  chamber  joined  to  the  nozzle  and  their  cal¬ 
culation  is  the  most  important  result  of  this  research  effort.*  The  cal¬ 
culations  have  been  performed  for  conical  nozzles  and  are  presented  in 
tabular  form.  Examples  are  presented  which  demonstrate  the  use  of  the 
tables  in  typical  problems.  Oscillatory  pressures  and  velocities  are 
also  calculated  in  a  limited  number  of  cases  in  order  to  provide  physical 
insight  to  the  oscillation.  An  asymptotic  analysis  has  been  performed 
which  is  valid  in  the  low  entrance-Mach-number  regime.  By  predicting  the 
admittance  coefficients  and  flow  properties  through  closed-form  solutions, 
the  asymptotic  analysis  is  an  asset  in  the  interpretation  of  results. 
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SOMNAIRE 


Un  traitement.  linearise  d’  un  e'coulement  oscillatoire  tridimensionnel 
dans  les  tuyeres  super critiques  a  ete  effectue  pour  les  deux  cas  ou 
1’ ecoulement  permanent  est  axisymetrique  ou  plan.  Dans  le  cas  axi- 
sym^trique  des  series  perturbatrices  ont  ete  utilisees  pour  l'^tude  des 
oscillations  non  lineaires.  Dans  les  analyses  effectuees  on  a  separe 
les  variables  en  reduisant  a  un  systeme  d’  equations  differentielles 
ordinaires  le  systeme  d’  equations  differentielles  partielles.  Les 
variables  qui  decrivent  les  dependences  transversales  des  propriety 
de  1* ecoulement  sont  regies  par  des  Equations  differentielles  connues 
(par  exemple,  on  obtient  1* equation  de  Bessel).  Par  contre,  les 
dependences  axiales  sont  decrites  par  des  equations  differentielles 
qu’ il  faut  resoudre  par  des  ooyens  numeriques.  Les  coefficients  d’  entree 
de  tuyere  sont  lies  aux  dependences  axiales  des  proprietes  de  1'  ecoulement. 
Certaines  techniques  ont  iti  appliquees  pour  reduire  l’ordre  des  Equations 
differentielles  en  vue  de  faciliter  le  calcul  des  coefficients  d’  entree. 

Ces  coefficients  d’  entree  seront  employes  pour  la  condition  limite 
appliquee  »  la  sortie  de  la  chaabre  relive  a  la  tuyere  et  leur  calcul 
repr^sente  le  resultat  le  plus  important  de  ces  efforts  de  recherches. 

Les  calculs  ont  ete  effsctues  pour  les  tuyeres  coniques  et  sont  presentes 
sous  forme  de  tableaux.  Des  exemples  sont  cites  pour  d£aontrer  1’  utilisa¬ 
tion  de  ces  tableaux  pour  la  solution  de  problemes  types.  Les  pressions 
et  les  vitesses  oscillatoires  sont  ^galement  calcul^es  dans  un  noabre 
limite  de  cas  pour  obtenir  une  connaissance  physique  de  1' oscillation, 
line  analyse  asymptotique  a  ete  realist  qui  est  valable  dans  les  conditions 
de  faible  noabre  de  Mach  d’  entree.  Elle  s’  avere  tres  avantageuse  pour 
1' interpretation  des  r^sultats  en  peraettant  de  prevoir  les  coefficients 
d’ entree  au  aoyen  de  solutions  de  forme  fera£e. 


POREVORD 


The  theory,  on  which  the  results  presented  in  this  monograph  are  based,  was  developed 
by  the  senior  author  over  a  decade  ago  and  presented  at  a  meeting  at  the  University 
of  Maryland.  At  that  time,  the  complete  formulation  of  the  linear  admittance  co¬ 
efficient  for  a  supercritical  nozzle  had  been  derived  in  the  general  ease  where  both 
vorticity  and  entropy  oscillations  exist  at  the  nozzle  entrance  In  the  following 
years,  a  small  number  of  calculations  were  made  with  the  purpose  of  providing  the 
necessary  data  for  the  determination  of  the  combustion  instability  limits  in  a  parti¬ 
cular  experimental  rocket11.  However,  the  publication  of  the  theory  was  postponed 
until  more  complete  calculations  and  interpretations  would  become  available. 

In  recent  years,  the  project  became  a  cooperative  venture  and,  primarily  through 
the  efforts  of  the  junior  author,  the  computer  calculations  and  other  related  analyses 
were  accomplished.  Meanwhile,  the  need  for  the  public  availability  has  become  pressing, 
due  to  the  problem  of  combustion  instability  in  rocket  motors.  One  of  the  leading 
rocket  manufacturers  in  the  United  States  has  found  it  necessary  to  establish  its  own 
computer  program  based  on  the  above-mentioned  theoretical  developments6.  Another 
scientist  has  decided  to  attack  the  problem  independently  in  a  somewhat-different 
tanner7.  Ve  welcomed  therefore  the  opportunity  (made  possible  by  AGARS)  to  publish 
the  theory,  some  of  its  extensions,  the  methods  of  calculation,  and  the  numerical 
results.  It  has  been  felt  that  they  would  fill  a  well-defined  void  in  the  existing 
technical  literature  and  play  a  useful  role. 


In  Part  I,  we  present  the  theoretical  background.  The  original  theory  of  the  senior 
author  is  contained  in  Sections  1  through  12.  Section  13  contains  the  discussion  of 
the  extension  of  the  theory  to  the  nonlinear  shockless  case  performed  by  B.T.Zinn, 
under  the  supervision  of  the  authors,  as  part  of  his  Ph.D.  thesis8.  Part  II  contains 
the  discussion  of  the  calculations,  related  analyses,  interpretation  of  the  results, 
and  examples  of  applications  and  is  mostly  due  to  the  junior  author.  Sections  14,  15, 
and  16  present  the  method  of  calculation  and  a  discussion  of  the  results.  Sections  17 
and  18  present  an  asyaptotic  theory  useful  for  the  purpose  of  interpretation.  Finally, 
Section  19  contains  a  few  sample  applications  of  the  results. 


Ve  wish  to  recognize  the  major  technical  support  provided  by  the  staff  of  the 
Guggenheim  Laboratories  Computing  Group  at  Princeton  University.  Also,  we  wish  to 
recognize  the  major  financial  support  provided  by  the  National  Aeronautics  and  Space 
Administration.  Additional  support  for  the  calculations  was  provided  by  the  William 
B.Reed,  Jr.  Fund  for  Engineering  Research  of  Princeton  University  and  (indirectly 
through  the  support  of  the  Princeton  University  Computing  Center)  by  the  National 
Science  Foundation. 
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BEHAVIOR  OF  SUPERCRITICAL  NOZZLES  UNDER 
THREE-DIMENSIONAL  OSCILLATORY  CONDITIONS 

Luigi  Crocco  and  William  A. Sirignano 


PART  I.  THEORY 


1.  INTRODUCTION 

Many  propulsive  devices  are  terminated  by  a  nozzle  through  shicb  ,te  propulsive 
gases  are  discharged.  Very  often  the  nozzle  operates  in  the  supercritical  range  and 
is  shaped  as  a  classical  Laval  nozzle,  converging  up  to  a  throat  (where,  in  steady- 
operation.  the  sonic  velocity  is  achieved)  and  diverging  thereafter. 

Unsteady  conditions  are  present  in  the  nozzle  when  the  operation  in  the  propulsive 
device  is  unsteady.  A  particular  type  of  unsteady  operation  which  has  great  iaportance 
in  practice  results  from  combustion  instability  in  the  propulsive  device.  In  this 
case  the  operation  is  oscillatory,  characterized  by  periodic  variations  of  the  flow 
parameters  both  in  the  combustion  chamber  and  in  the  nozzle  which  follows. 

The  study  of  combustion  instability  is  of  great  importance  for  the  safety  of  opera¬ 
tion  of  combustion  devices,  and  has  been  already  the  subject  of  a  great  deal  of  theore¬ 
tical  and  experimental  research.  It  is  important  in  these  studies  to  hnow  the  behavior 
of  the  nozzle  under  oscillatory  conditions.  In  particular,  it  i3  necessary  to  find 
out  hn«  a  wave  generated  in  the  combustion  chamber  is  partially  reflected  and  partially 
transmitted  at  the  entrance  of  the  nozzle.  Mathematically,  this  is  equivalent  to 
saying  that  it  is  necessary  to  know  the  hwmdary  conditions  created  by  the  nozzle  to 
the  oscillatory  flow  present  in  the  combustion  chamber.  For  instance,  on  any  solid 
wall  of  the  chamber,  the  boundary  condition  is  that  the  velocity  component  of  the 
gases  normal  to  the  wall  must  be  zero  at  every  instant,  and  so  must  be.  therefore,  the 
corresponding  velocity  perturbation,  no  natter  what  are  the  perturbations  of  the 
tangential  velocities,  pressure,  entropy,  etc.  Obviously  if  the  entrance  of  the  nozzle 
is  considered  to  represent  one  of  the  boundaries  of  the  combustion  chamber,  the  cor¬ 
responding  boundary  conditions  are  more  involved.  No  single  perturbation  can  be 
assumed  to  vanish,  and  the  boundary  condition  may  be  expected  to  be  expressed  by  a 
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relation  between  the  various  perturbations  which  shall  be  called  the  admittance 
condition  .  In  particular,  if  the  perturbations  are  assumed  to  be  of  sufficiently 
small  amplitude,  so  that  the  problem  can  be  linearized,  the  aforesaid  relation  express¬ 
ing  the  boundary  condition  must  be  linear.  Observe  that  a  linearized  treatment  can 
be  only  applied  to  the  study  of  incipient  instability,  that  is  to  the  study  of  the 
combustion  stability  limits  of  a  given  propulsive  device.  When  the  instability  is 
fully  developed,  nonlinear  effects  become  essential  and  must  be  taken  into  consideration. 

If  the  unperturbed  flow  in  the  nozzle  can  be  assumed  to  be  one-dimensional,  a  parti¬ 
cularly  simple  case  is  obtained  when  also  the  perturbations  are  taken  one-dimensional, 
that  is  the  perturbations  are  assumed  to  be  uniform,  at  each  instant,  on  each  section 
of  the  nozzle.  Obviously  in  this  case  the  conditions  in  the  combustion  chamber  are 
also  one-dimensional,  and  any  wave  present  in  the  chamber  or  in  the  nozzle  is  of  an 
“axial"  type.  The  behavior  of  the  nozzle  in  the  presence  of  axial  waves  has  been 
analyzed  by  Tsien1  in  a  few  simple  cases,  and  by  Crocco2’3  for  the  most  general  type 
of  linear  axial  oscillation,  in  particular  when  the  entropy  is  not  constant.  The 
variability  of  the  entropy  is  in  fact  an  unavoidable  consequence  of  the  combustion 
taking  place  under  oscillatory  conditions,  and  may  product  interesting  effects  on  the 
combustion  instability3.  Exy.i iments  carried  out  some  time  ago  in  nearly  isentropic 
conditions  show  satisfactory  agreement  with  the  theoretical  predictions4. 

However,  the  axial  type  of  oscillations  is  only  a  particular  case,  and  in  actual 
conditions  transverse  perturbations  are  often  present,  this  being  particularly  true 
in  large  propulsive  devices.  The  problem  is  now  complicated  by  the  incre'  'ed  number 
of  degrees  of  freedom,  but  under  the  same  assumptions  that  the  perturt  -tious  are 
small  and  that  the  unperturbed  flow  in  the  nozzle  is  one-dimensional  (hence,  irrota- 
tional),  isoenergetic  and  isentropic,  it  is  amenable  to  a  relatively  simple  analytical 
treatment.  It  is  the  purpose  of  this  monograph  to  show  how,  in  this  simpler  case, 
the  admittance  condition  at  the  nozzle  entrance  can  be  expressed  and  to  discuss  the 
corresponding  numerical  results.  We  shall  treat  both  the  case  of  axisymmetric  nozzle 
(the  most  common  in  practical  devices)  and  that  of  two-dimensional  nozzle.  We  shall 
also  briefly  discuss  the  nonlinear  treatment  of  the  nozzle  admittance  problem  in  the 
absence  of  shock  waves. 

Observe  that  the  results  of  the  present  study  are  applicable  to  any  type  of  device, 
propulsive  or  not,  involving  combustion  processes  or  not. 

*  The  following  considerations  may  help  in  understanding  the  nature  of  the  admittance  condition. 

It  is  clear  that,  if  the  flow  in  the  nozzle  is  supercritical,  for  sufficiently  small  oscilla¬ 
tions  the  supersonic  portion  of  the  nozzle  has  no  effect  on  the  chamber  conditions,  because 
downstream  of  the  throat  the  oscillations,  no  matter  how  distorted,  must  always  propagate 
downstream  and  cannot  interfere  wlch  the  upstream  flow.  Hen re,  the  logical  choice  for  the 
surface  on  which  boundary  conditions  must  be  prescribed  would  be  the  surface  "here  "ie  sonic 
velocity  is  achieved,  or,  for  small  oscillations  around  an  approximately  one-uimensional  flow, 
the  throet  Itself.  It  has  been  shown2-3  that  the  proper  boundary  condition  at  the  throat  is 
that  the  solution  remains  regular  here  (where,  Indeed,  a  singularity  tends  to  result  from  the 
Inability  of  the  disturbances  to  propagate  upstream  from  the  supersonic  into  the  subsonic 
region).  In  practice,  however,  it  is  useful  to  divide  the  whole  of  the  chamber  plus  the 
nozzle  in  two  parts:  the  combustion  chamber  extending  down  to  the  nozzle  entrance  where  the 
processes  of  combustion  are  taking  place  but  the  mean  flow  Mach  nuaber  is  relatively  low; 
and  the  nozzle  where  no  combustion  is  assumed  to  take  place  but  the  mean  Mach  number  grows  up 
to  unity.  The  result  of  this  subdivision  is  to  move  the  boundary  of  the  combustion  chamber 
from  the  throat  up  to  the  nozzle  entrance,  where  the  appropriate  boundary  conditions  can  be 
obtained  by  studying  the  oscillatory  behavior  of  the  nozzle  per  te  and  obtaining  the  proper 
relation  between  the  perturbations  at  the  nozzle  entrance  (the  admittance  condition)  from  the 
condition  of  non-singularity  at  the  throat. 
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2.  THE  EQUATIONS 

Using  stars  to  denote  dimensional  quantities  and  operators,  the  equations  of  motion 
for  an  invlscid,  non-heat-conducting  gas  are  the  following. 

Continuity: 

BP* 

— -  +  V*.  (pV)  =  o  . 

o  t 


Momentum: 

~  +  ~  V(q*2)  +  (V*  x  q*)  x  q*  =  -  -I  W  . 
at*  2  p* 

where  t*  represents  the  time,  p*  and  p*  the  density  and  pressure,  and  q*  the 
velocity  vector. 

Wien  viscosity  and  heat  conductivity  are  disregarded  the  energy  equation  in  its 
simplest  form  expresses  the  constancy  of  the  entropy  for  a  fluid  particle  after  it 
enters  the  nozzle: 

~  +  q*.V  V  =  0  , 
dt 


where 


„★  - 


cD  (-  log  p*  -  log  p*)  +  constant 


(l) 


represents  the  entropy.  The  specific  heat  cp  is  assumed  to  be  constant.  With  this 
particular  form  of  the  energy  equation  it  is  not  necessary  to  introduce  the  equation 
of  state,  which  is  implicitly  taken  into  account  in  Equation  (1).  and  the  four  equations 
just  written  are  complete  in  the  unknowns  q*  ,  p*  ,  p*  .  s*  . 


These  equations  can  be  nondimensionalized  using  appropriate  reference  values. 
Assuming  the  gas  entering  the  nozzle  in  the  unperturbed  flow  to  be  isoenergetic  and 
isentropic  (as  well  as  irrotational),  and  hence  to  stay  such  in  the  following  expansion 
through  the  nozzle,  the  corresponding  stagnation  quantities  remain  constant  throughout 
the  unperturbed  flow  and  hence  are  suitable  reference  quantities.  Hence  we  define 


q 


(2) 


where  c  denotes  the  sonic  velocity,  the  superscript  0  indicates  stagnation  values 
and  the  superposed  bar  unperturbed  (steady)  values.  Observe  that 


*  -  **.  -w-  u-*- *«*  *  m- •  nw*'  safe-  «mm»i 
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The  lengths  can  be  nondimenslonalized  using  a  suitable  characteristic  length  L*  to 
be  further  defined,  and  a  nondimensional  time  is  immediately  obtained  as 


t 


(4) 


Transforming  also  the  operators  to  nondimensional  coordinates  we  obtain  the  equations 
of  motion  in  the  form 


}P 

Bt 


+  V.(pq)  =  0 


— —  +  -  V(q2)  +  (Vx  q)  x  q  = - Vp 

Bt  2  yp 


Bs 

— —  +  q.Vs  =  0  , 
Bt 


with 


s*  1 


s  =  —  =  -  log.  p  -  log.  p  t  constant 


representing  the  nondimensional  entropy.  Equation  (3)  is  replaced  by 


c°* 


* 


c  representing  the  nondimensional  sonic  velocity. 


(5) 

(6) 

(7) 

(8) 

(9) 


3.  LINEARIZATION  OF  THE  EQUATIONS 

We  now  introduce  the  ordinary  assumption  that  the  (unsteady)  perturbations  around 
the  (steady)  unperturbed  quantities  are  of  such  small  amplitude  that  only  linear  terms 
in  the  perturbations  need  be  considered.  Hence  we  introduce 

q  =  q  +  q'  ;  p  =  p  +  p'  ;  p  -  p  +  p'  ;  s  =  5  +  s'  ( 10) 

in  Equations  (5)  to  (8)  and,  after  neglecting  all  terms  of  order  higher  than  first  in 
the  primed  quantities,  we  separate  each  equation  into  its  unperturbed  portion  (containing 
only  the  barred  quantities)  from  the  perturbation  portion  (linear  in  the  primed 
quantities).  Clearly  each  portion  of  the  equations  must  be  satisfied  separately.  The 
equations  for  the  unperturbed  flow  are 


grad  —  +  (V  x  q)  x  q  =  - — grad  p 

2  yj> 


div  (/3q)  =  0  ; 


l 
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q.grad  s  =  0  ;  s  =  -  log  p  -  log  p  +  constant  . 


These  equations  can  be  replaced  by  the  following  simpler  ones,  obtained  from  (11)  in 
a  standard  fashion  when  the  flow  is  irrotational, 

div(£q)  =  0;  D  =  Py  ;  4=1  - - “  Q2  •  (12) 

P  2 

The  second  expresses  the  constancy  of  the  entropy  s  ,  the  third  that  of  the  stagnation 
temperature.  Comparing  this  last  equation  to  (9)  we  obtain  the  unperturbed  sonic 
velocity 

c7  =  1  -  — - -  q7  =  p<y-1'>  ,  (13) 

2 

the  last  term,  resulting  from  (12),  providing  the  unperturbed  density.  Ihen  the  nozzle 
is  axi symmetric  or  two-dimensional  the  first  (12)  can  be  used  to  define  a  stream 
function.  We  may  write  indeed 

r£q  =  e0  x  V</»  (14) 

in  the  axisymmetric  case  (r  representing  the  nondimensional  distance  from  the  axis  of 
symmetry  and  e0  the  unit  vector  in  the  tangential  direction)  and 

bpq  =  ey  x  Vi p  (15) 

in  the  two-dimensional  case  (b  representing  the  nondimensional  width  of  the  nozzle, 
which  must  be  assumed  to  be  finite  if  three-dimensional  oscillations  are  to  be  con¬ 
sidered,  and  ey  the  spanwise  unit  vector.) 

In  the  present  assumption  of  irrotationality  of  the  unperturbed  flow,  a  potential 
function  can  also  be  defined  as 

Q  =  V0  .  (16) 

The  stream  and  potential  functi  /us  introduced  are  nondimensional. 

In  what  follows  the  unpe-.  turbed  flow  will  be  assumed  to  be  in  the  meridional  plane 
(axisymmetric  case)  or  purely  two-dimensional  (two-dimensional  nozzle). 

The  perturbation  equations  can  be  written  as 


+  V.  (qp'  +pq')  =  0 


Bq'  Ip 1 

— —  +  7(q.q')  +  (V  x  q')  x  q  - - ;Vp'  +— — Vp 

ot  yp  ypi 
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+  q.Vs'  =  0 


n 1  p1 

s’  =  —  ~  -  (20) 

yv  p 

Equations  (18)  to  (20)  represent  the  system  of  equations  to  be  solved.  It  is  linear 
in  the  perturbations,  and  has  coefficients  depending  on  the  solution  of  the  unperturbed 
Equations  (12). 

It  is  useful  to  rewrite  the  second  member  of  (18)  in  a  different  form,  using  the 
following  transformations: 

1  p'  (d'\  p'  p' 

— -Vv'  +-WVp  =  -V  -hr  --X.VP+-WVF 

yp  yp  \ypj  yp  yp 


Vp/p'  p' 


yp  yp  \yv  ~P. 


-  -  vf-4)  +  Kv?)s'  , 

\ypj 


where  use  has  been  made  of  the  second  (12),  of  the  second  (11)  and  of  (20). 


4.  CHOICE  OF  THE  INDEPENDENT  VARIABLES:  A’lSYHHETRIC  NOZZLE 

Abandoning  the  vectorial  representation,  it  is  useful  to  choose  the  Independent 
spatial  variables  in  a  way  appropriate  to  the  introduction  of  thb  boundary  conditions 
at  the  nozzle  walls.  In  the  case  of  axial  symmetry  a  suitable  choice  is  to  take  the 
steady-state  potential  function  to  replace  the  axial  variable,  and  the  steady-state 
stream  function  to  replace  the  radial  variable.  Indicating  with  5s  and  8n  ele¬ 
mentary  (nondimenslonal)  lengths  in  the  direction  of  the  unperturbed  streamlines  and 
of  their  normal  in  the  meridional  plane  of  Figure  1,  Equations  (14)  and  (16)  can  be 
written 

d  <p  __  d'p 

5=7^;  rAa  =  (22) 

os  on 


Hence  the  *7-  re  of  the  elementary  length  dl  can  be  written,  in  terns  of  d</>  .  d</»  , 
and  dd  (6  representing  the  azimuthal  variable)  as 

dl2  =  a |<tf2  +  h|d 4?  +  h|d<92  =  4fdtf>2  +  d'p2  +  r2  dd*  , 


from  which  we  get 


1  1 
h,*  -  _  >  hj,  - 
9  q  ^  rpq 


he  =  r 


(23) 
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These  quantities  are  used  to  calculate  the  divergences,  gradients,  and  rotors  appearing 
in  the  equations.  If  f  is  any  scalar  and  F  =  F^e^  +  F^e^  +  Feee  is  any  vector, 
e^  ,  e^  .  ed  being  unit  vectors  normal  to  the  surfaces  0  =  constant  ,  0  =  constant  , 
6  =  constant  ,  we  have 


__  /  a  v  x 

Vf  =  — — 


_ii£ 


__  Bf 


e*+(rJ)e* 


=  °^}  (h*W  +Te  (W*>_ 

=  3/vr 

B0\^<y  30  \  q  J  3#  \rpqy 


K*4> 

Vs* 

1 

3 

3 

3 

h<t>b^6 

30 

30 

3? 

W 

k|/,fyr 

b6?6 

=  Pq 


3  ,  D  \  B  Q  3  {V*\  ^  ,  „  . 

ff ^  ^  e* + ?  |p  (ij  -  *  (rp^_ 

B0  \rP<y  30  \  q/ 


Using  these  relations,  the  conservation  equations  can  be  written  explicitly,  noticing 
that  q  =  qe^,  and  defining  the  components  of  the  velocity  perturbation 

q'  =  u'e^  +  v'fy  +  w 'e0  .  (27) 

We  obtain: 

Continuity  (Equation  (17)  divided  by  p): 


—  (^\  +  (—  +  —  ]  +  (r*  2  v'  \  +  _L 

Bt  \p)  11  30  Ip  q/  Q  30  V  ^  r pq/  rz  'd& 


Momentum  Equation  (18),  taking  into  account  (21): 


^-component  (divided  by  q): 

3  /u'\  +  3  /-j  u,N\  3  /p'\  _  1  3q2  t 

3t  \q/  30  y  q )  30\yp/  2  30 

0-component  (divided  by  r/33): 

3  /  v'  \  —  3  /v'  \  3q2  u'  3  /p'\  lBq2  f 

dt\rpqj  Q  30\rpq/  30  q  30  \yp/  2  30 


(29) 


(30) 


0-component  (multiplied  by  r)  taking  into  account  that,  in  view  of  its  axial 
symmetry,  the  unperturbed  flow  does  not  depend  on  0  : 


U  — T  u  o  /  p 

—  (rw')  +?-(«')  +—  -= 

3t  30  30  Wp 


Entropy,  Equation  (19): 


3s'  —y3s' 

-  +  q‘ - 

-3t  30 


0  . 


In  view  of  (9),  Equation  (20)  can  be  rewritten  in  the  form 


c7  TP  p 


(32) 


(33) 


The  preceding  equations  have  been  purposely  written  in  a  way  that  shows  that,  instead 
of  J  ust  the  perturbations,  it  is  convenient  to  choose  the  dependent  variables  to  be 
the  combinations  p'/p  ,  u'/q  ,  v'/rpq  ,  rw'  ,  p'/yp  ,  and  s'  . 

We  see  at  once  that,  from  the  preceding  six  equations  in  these  six  dependent  variables, 
one  can  obtain  a  system  of  four  equations  in  the  first  four  of  the  dependent  variables 
just  listed.  In  fact,  s'  can  be  obtained  from  Equation  (32)  independently  of  the 
other  equations  once  the  corresponding  s '-distribution  is  prescribed  at  the  nozzle 
entrance  as  a  function  of  time,  after  which  p'/yp  can  be  obtained  from  (33)  in  terms 
of  p'/p  .  However,  the  resulting  system  is  still  way  too  complicated  to  be  amenable 
to  solution.  A  major  step  in  the  way  to  a  solution  would  be  accomplished  if  the  vari¬ 
ables  were  separable,  which  they  are  not  in  the  system  as  it  is. 


S.  CHOICE  OF  THE  INDEPENDENT  VARIABLES:  BIDIMENSIONAL  NOZZLE 

If  the  nozzle,  and  the  corresponding  unperturbed  flow,  are  two-dimensional,  the 
independent  variables  are  chosen  in  a  similar  way.  According  to  Equations  (15)  and 
(16),  Equations  (22)  can  be  replaced  by 
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</»- component  (divided  by  bpq): 

B  /  v'  \  -5  3  /  v'  \  B  q2  u'  B  / p'^  1  Bq2 

Bt  \bp  q/  Q  B'/'Vbpq/  By  q  a4‘\7Pj  2  3^ 


y-component,  taking  into  account  that  the  unperturbed  flow  does  not 

Bw'  -5  Bw'  B  / p'\ 

-  +  q2 -  + — (—  =  0  . 

Bt  By  By\yp/ 


depend  on  y 

(41) 


Entropy: 


Bs'  ,  Bs' 

s'  -  -L£l 

Bt  Q  By 

■  c7  yp 

(42) 


These  last  two  equations  are  the  same  as  for  the  axisymeetric  nozzle,  and  can  be  used 
to  reduce  to  four  the  number  of  dependent  variables.  However  the  same  difficulty  as 
before  is  found,  due  to  the  complexity  of  the  equations  which  prevents  the  separation 
of  the  variables. 


6.  SEPARATION  OF  THE  VARIABLES  FOR  ONE-DIMENSIONAL 
UNPERTURBED  FLOH:  AXISYMMETRIC  NOZZLE 

The  chief  obstacle  that  prevents  the  separation  of  the  variables  in  Equations  (28) 
to  (33)  is  the  fact  that  q  and  p  depend  on  both  4>  and  ^  in  a  way  determined 
by  the  solution  of  (12).  An  additional  obstacle  is  due  to  the  presence  of  the  factors 
r2  and  r'2  appearing  in  the  two  last  terms  of  Equation  (28).  These  factors  should 
be  expressed  in  terms  of  the  independent  variables  4>  ,  $  in  a  way  depending  again 
on  the  solution  of  (12). 

Botn  obstacles  are  removed  if  the  unperturbed  flow  is  one-dimensional.  This  means 
that  the  dependence  of  q  and  p  on  y  can  be  practically  disregarded,  so  that  they 
can  be  considered  practically  uniform  on  each  surface  <p  -  constant  .  It  means  also 
that  the  angle  of  obliquity  of  the  streamlines  with  respect  to  the  axis  of  symmetry 
is  sufficiently  small  so  that  its  cosine  is  practically  1  and  the  element  of  normal 
Sn  along  the  surface  4>  =  constant  can  be  identified  with  dr  .  Hence  Equation  (22) 
can  be  integrated,  providing 

*  •= 

or 


which  provides  a  simple  expression  for 
advantage  of  dropping  the  dependence  of 
vanish  identically. 


p(0)q(0) 


2*p 
P  q 


('.3) 


r2  ,  to  be  introduced  in  (28).  An  additional 
q  on  is  ?hat  two  terms  of  Equation  (30) 


Without  writing  down  explicitly  the  equations  obtained  assuming  q  =  q(<£)  ,  p  -  p(<t>) 
and  introducing  (43),  we  observe  that  the  variables  are  now  separable.  In  fact,  if 
we  assume  a  harmonic  time  dependence,  expressed  in  the  complex  form,  and  express  the 
dependent  variables  as  follows*: 

—  =  R(0W>Q(S)eIiJt 

P 


=  v(<pm^)melut 


—  =  V^f'^We1*1 

rpq 


=  W(0W0)0'(6?)e 


i /as*,  i^t 


p  cmmmt 


s'  =  s(0m</oe(0)ela,t  ,J 

with  ¥'(</»)  =  d'Ji/d'/'  ,  8'  {6)  =  d8/d#  .  and  u>  representing  the  nondiaensional  angular 
frequency,  related  to  the  diaensional  co*  by 


CO  =  ~co 

„o* 


The  equations,  divided  by  ?8  exp  (tcot)  ,  take  the  following  form: 


-j  ,  -  ,  V'  ?\  PPW  0" 

+  qR  Vu +  2^v  (*  7  +  +-^"e  = 


d  .  ,1  dq 

iwu  + — (q2U)  +  P' - -S  =  0 

d0  2  d<£ 


(i«V  +  q*  V'  +  P)  —  =  0 


(iwW  +  qV  +  P)  —  =  0 
o 


#  The  various  functions  appearing  in  (44)  may  be  complex  functions  of  the  corresponding  variable. 
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iwS  +  qV  =  0  (50) 

S  =  4-P  -  R  ,  (51) 

cz 

where,  again,  primes  represent  differentiations  with  respect  to  the  corresponding 
independent  variable. 

It  is  interesting  to  express  the  coaponents  of  the  flow  vorticity  using  (26)  and 
(44).  The  result  is 


V  x  q  =  pq(W  -  V)?'0'ei"te(i)  +  —  (U  -  W')90,e1"t«^  +  rpq2(V'  -  U)f'8elc,,te0  .  (52) 

We  see  that  the  variables  are  separated  in  Equations  (47)  to  (51),  which  are  made  up 
of  factors  depending  on  a  single  independent  variable.  Thus  the  only  requirement  for 
complete  separation  is  that  the  variables  and  6  disappear  from  (46).  For  8 
to  disappear.  9"/8  must  be  a  constant;  more  precisely,  in  view  of  the  necessary 
periodicity  of  9  ,  we  must  have 


9" 

T 


-  v 


(53) 


v  representing  an  integer. 

Introducing  (53)  in  Equation  (46),  it  appears  that  the  dependence  on  > p  disappears 
only  if  V  =  f  and  ’?(</')  satisfies  the  equation 

yn  ijj/  v2 

*P  ~zr  +  ■nr - r  =  constant  .  (54) 

®  t  4</> 

with  the  value  of  the  constant  to  be  determined  from  the  boundary  conditions. 

This  result  shows  that  the  requirement  for  the  separation  of  the  variables  introduces 
a  certain  restriction  on  the  solution  of  our  equations.  In  fact,'  the  condition 

V  =  W  (55) 

is  equivalent,  as  (52)  shows,  to  the  condition  that  the  vorticity  component  of  the 
perturbed  flow  along  the  streamlines  must  be  identically  zero.  The  resulting  loss  of 
generality  does  not  appear  too  important,  since  in  most  practical  cases  no  mechanism 
is  present  to  generate  axial  vorticity  in  combustion  chambers. 

The  general  solution  of  (53)  is* 

9  =  cos  v(8  -  a)  ;  9  =  ^M^-a)  t  (56) 

where  a  represents  a  constant,  and  an  unessential  multiplicative  constant  has  been 

taken  to  be  unity. _ 

*  The  following  discussion  concerning  the  6  and  ?  functions  proceeds  along  the  well  known 
lines  used  for  acoustic  oscillations  in  cylinders^.  Ve  repeat  It  here  for  the  purpose  of 
completeness. 
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The  first  (56)  leads  to  a  standing  mode  of  oscillation,  as  clearly  seen  froc  the 
fact  that  the  aaplitude  of  each  of  the  perturbations  defined  hy  (44)  vanishes  r-n  v 
fixed  diametric  il  planes  of  the  nozzle.  The  second  (56)  leads  to  a  rotating  BGde  of 
oscillations,  where  the  aaplitude  of  the  perturbations  reaains  constant  for  fixed 
values  of  ±  vd  +  cot  .  The  angular  speed  of  the  rotating  Bodes  is  d#/dt  =  .  In 

both  cases  v  represents  the  nuaber  of  the  diaaetral  nodal  planes  (standing  or  rotating) 
for  the  node  under  consideration. 

Next,  consider  (54).  Taking  (sW*  as  independent  variable,  this  is  transforaed 
into  a  Bessel  equation  of  order  v  .  The  general  solution,  if  the  constant  appearing 
on  the  right-hand  side  is  negative*,  and  set  equal  to  -  k2/4  ,  is  given  by 

*  -  C3Jv(k[$*>  +  C,Yv(k[</0*)  . 

where  J  and  Y„  are  Bessel  functions  of  the  first  and  second  kind,  respectively. 

Since  Yv  does  not  reteain  finite  as  <p  approaches  zero.  Ca  aust  vanish  identically 
if  the  perturbations  are  to  be  finite  on  the  nozzle  axis.  Hence,  setting  the  un¬ 
essential  Bultiplicative  constant  C3  to  be  unity,  we  have 

5(0  =  J„(k[#)  .  (57) 

The  constant  lr  is  detersined  by  considering  the  boundary  condition  v'  =  0  at 
the  nozzle  wall,  which  is  a  =  constant  surface.  Calling  >/',  the  value  of  the 
streaa  function  at  the  wall*,  and  recalling  (43),  the  condition  v'  =  0  at  the  wall 
is  satisfied  if 

J.'.(k£<pj«)  =  0  .  (58) 


The  priae  here  indicates  differentiation  with  respect  to  the  argument.  Hence,  if 
svb  is  the  htJ>  zero  of  J.'(x)  .  in  order  to  satisfy  (58)  we  oust  have  k  =  svb/(</'f)* 
and  (57)  becomes**: 


(59) 


the  last  tara  being  obtained  froa  (43).  Observe  that  the  last  for*  of  ¥(i/0  and  the 
forms  (56)  of  @(£)  are  exactly  the  sane  as  for  the  acoustic  oscillations  of  a  gas 
in  a  cylindrical  chamber,  the  only  difference  being  that  in  (59)  rw  is  the  variable 
radius,  rM(<p)  ,  of  the  nozzle  sections,  instead  of  a  constant.  Hence  within  the 
validity  of  the  one-dimensional  assumption  for  the  unperturbed  flow,  the  perturbations 
are  distributed  on  each  section  of  the  nozzle  in  the  same  way  as  they  are  in  the  acoustic 
oscillations  in  a  cylindrical  chamber.  Corresponding  to  the  lowest  values  of  u  and 
h  we  have  for  s„h  the  values  given  in  Table  I  on  p.84  (Ref.  5)**. 

*  Positive  values  of  the  constant,  leading  to  the  Bessel  fractions  of  iaaginary  argument,  would 
not  allow  the  boundary  condition  at  the  nozzle  wall  to  be  satisfied. 

t  The  value  of  ^  Is  laeodiaicly  obtained  froa  (43)  evaluated  at  the  throat  or  at  the  nozzle 
entrance,  once  the  reference  length  L*  has  been  chosen.  Appropriate  values  for  L*  are 
either  tbo  radius  of  the  throat  section  or  the  radius  of  the  entrance  section  of  the  nozzle. 


**  The  index  h  has  a  aesning  siwiiar  to  that  of  v  ;  b  -  1  represents  the  nuaber  of  nodal 
circles  for  the  pressure  perturbations,  that  is.  the  nuaber  of  circles  on  each  section  on 
which  the  pressure  does  not  oscillate. 


H 


Ihs  arguaent  of  the  Bessel  function  in  Reference  5  is  "s^r/r,  so  that  tbetr  eigenvalues 
aust  be  wultiplted  by  T'  in  order  to  obtain  the  values  In  Table  I. 


or  j aan »  yo-wiji  mm  hn— 


■rBwetfac 
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Inserting  in  Equation  (46)  the  relations  (55),  (53),  and  (54)  with  the  value  of  the 
constant  in  the  last  relation  given  by  -  k2/4  =  -  s2h/4</'w  ,  we  obtain 

_  _  g2 

icoR  +  qV  +  q2u'  -  -^pqV  =  0  .  (60) 


The  separation  of  the  variables  is  now  complete. 


7.  SEPARATION  OF  THE  VARIABLES  FOR  ONE- DIMENSIONAL 
UNPERTURBED  FLOW:  BIDIMENSIONAL  NOZZLE 

When  q  and  p  are  taken  to  be  functions  only  of  4>  .  Equations  (38)  to  (42) 
become  separable  too.  We  take 

~  =  R1(0)'?,(>/')Y(y)eia’t 


—  =  U1(0)?1(<//)Y(y)eI 


4=  =  V1(0)'P1'(^)Y(y)e1"t 
>pq  1 

w'  =  W1(0)f1(</»)Y,(y)ei&,t 

~  =  P1(<A)'i'1('/')Y(;')ela’1' 
yp  1  X 

s'  =  S1(0)?,1(v)Y(y)ei"t  , 


and  replace  these  expressions  in  the  equations.  Dividing  all  equations  by 
'PjO/OzU)  exp  (icixt )  ,  we  obtain 

icoRl  +  ?R[  +  +  bVq2V1  +  W,  Y  =  0 


-s  dq2  1  dq2 

i,  +  qV.  +  — o.+p - -  s.  = 

1  1  d0  1  1  2  d0  1 


(1«V1  +  ?vi  +  Pj)  =  0 


(icoWj  +  q*Wj  +  Pt)  —  =  0 


iwSj  +  d2S[  =  0 


s!  =  ffPi  "Ri  • 
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(66) 

(67) 


The  flew  vorticity  is  given  by  (37)  after  substitution  of  (61): 

V  x  <i  =  bpq(Wj  -  V^yV^  +  q(Uj  -  W')¥1Y/e1*tfy  + 

+  bpq^(Vj  -  U^Ye1"** 


(68) 


We  see  that  in  Equations  (63)  to  (67)  the  variables  are  separated, 
separation  may  be  complete  we  must  have 


¥ 


-  =  C, 


Y" 

—  =  Cy  , 
Y  y 


In  order  that  the 


(69) 


Cy  and  Cy  representing  constants  to  be  determined  through  the  boundary  conditions. 

We  notice  that,  contrary  to  the  axisymmetric  nozzle,  in  the  two-dimensional  case 
no  restriction  is  placed  on  the  solution  by  the  requirement  of  variable  separation, 
and  that,  according  to  (68),  the  vorticity  of  the  flow  may  have  non-vanishing  components 
in  any  direction.  The  values  of  C ^  and  Cy  are  determined  by  the  boundary  conditions 
of  vanishing  normal  velocity  at  the  nozzle  walls,  that  is  v'  =  0  at  0=0  and 
0  =  0W  ,  and  w'  =  0  at  y  =  0  and  y  =  b  *.  We  observe  that,  since  in  the  assump¬ 
tion  of  small  obliquity  8n  may  be  identified  with  dx  in  the  second  Equation  (34), 
x  representing  the  nondimensional  coordinate  along  the  nozzle  height,  we  obtain 


0  =  b/5qx  ,  0  =  bjoqa  , 


(70) 


where  a  =  a (0)  is  the  nondimensional  height  of  the  nozzle  section  under  consideration. 
The  value  of  0W  (corresponding,  as  for  the  axisymmetric  nozzle,  to  the  total  flux 


of  mass  through  the  nozzle)  can  be  determined  once  the  reference  length 
chosen.  According  to  (61)  the  boundary  conditions  can  be  written 

?J(0)  =  ?|(0W)  =  0  ;  Y'(0)  =  Y'(b)  =  0  . 

The  proper  solutions  of  (69)  satisfying  (71)  are 


(0)  =  cos 


Y(y)  =  cos  ( N77  - 


has  been 


(71) 


(72) 


*  For  the  sake  of  syaaetry,  the  ’‘axis”  of  the  two-dloensional  nozzle  should  be  taken  as  the 
line  0  =  0  ,  y  =  0  .  It  is  acre  ci-uvenient,  however,  to  take  0=0  on  one  of  the  curved 
walls  of  the  nozzle,  and  y  =  0  on  one  of  the  plane  walls.  The  other  curved  wall  is  then  a 
0  =  constant  =  0W  surface,  and  the  other  plane  wall  is  the  y  =  b  surface. 
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M  and  N  representing  integral  numbers  defining  the  mode  of  transverse  oscillation. 
Their  meaning  is  immediately  found  to  be  the  number  of  nodal  surfaces  for  the  pressure 
perturbation  in  the  corresponding  direction.  Equations  (72)  are  the  same  as  for  the 
acoustic  oscillations  of  a  gas  in  a  chamber  of  rectangular  cross-section,  so  that  again 
we  find  that  under  the  assumption  of  one-dimensional  unperturbed  flow  the  perturbations 
are  distributed  in  a  similar  fashion  on  each  cross-section  of  the  actual  duct  as  they 
are  in  the  acoustic  oscillations  in  a  rectangular  chamber. 

The  comparison  of  (72)  and  (69)  gives 

M  V 

\  =  C*  =  "  : 

Replacing  these  values  in  (62),  and  writing 

M7Tb 


we  obtain 

iwRi  +  q2R'  +  q2Uj  -  S2pVVj  -  cr2Wj  =  0  .  (74) 

The  separation  of  the  variables  is  now  complete. 


8.  REDUCTION  OF  THE  SYSTEM:  AXISYXMETRIC  NOZZLE 

We  first  observe  that  the  case  v  =  0  ,  h  =  1  ,  resulting  in  s,^  =  0  corresponds 
to  purely  axial  oscillations,  since  0  and  ?  become  constants.  In  this  case  (48) 
and  (49)  are  automatically  satisfied,  and  Equations  (46),  (47),  (50),  and  (51)  contain 
only  the  variables  R  ,  U  ,  P  ,  and  S  .  The  solution  for  this  case  has  been  discussed 
in  References  2  and  3,  but  it  will  again  be  included  as  a  particular  case  in  the  follow¬ 
ing  treatment. 

Similarly,  the  cases  in  which  either  v  =  0  or  h  =  1  (purely  radial  or  purely 
tangential  oscillations)  can  be  included  in  the  following  treatment  as  particular 
cases.  We  shall,  therefore,  discuss  the  case  where  v  and  h  are  arbitrary  integers. 

In  this  case,  in  order  to  satisfy  (48)  and  (49),  the  corresponding  expressions  in 
parentheses  must  vanish.  Subtracting  one  from  the  other  we  eliminate  P  and  obtain 


Y" 

—  =  Cv 
Y  y 


nV 


for  brevity 

N77 

°n  =  * 


(73) 


— t  d 

iw(W  -  V)  +  q2  — -  (W  -  V)  =  0  . 
A<p 

The  general  solution  of  this  can  be  written  in  the  form 


(75) 


=  c0f0 


W  -  V 


(76) 


where,  for  brevity,  *e  have  written 
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with 


t<t> 


i2(<P)  =  fo(0) 


7  «>')f  „(<*')  0 


(83) 


th 


Using  (78),  Equation  (80)  can  be  rewritten  as 

icof  +  q2$'  +  c2R  =  a(ll  -  c2f0)  . 

Equation  (60)  can  be  written  in  the  sane  dependent  variables  making  use  of  (82),  as 

2  2 

IcjR  +  qfy  +  q2$"  -  ^-pq$  “  --^p-p qfcrf2  +  C.f0)  . 

2  <pw  2</'w  1  10 


Elimination  of  R  from  these  two  equations  and  use  of  (13)  produce  the  following 
equation 


7(7  -  7)$"  -  7 


1  d7 

(7  d0 


+  2iw  $'  +  \or  - 


2  7  -  1  io>q2  dq2  s2h - ? \  * 


XT? 


ico 


d<£  \c 


2'A» 


-  c,  -^pqC2fn  .  (84) 


2*. 


This  is  the  fundamental  equation  to  which  our  system  can  be  reduced;  its  solution  can 
only  be  obtained  by  numerical  integration*,  using  the  appropriate  boundary  conditions. 

The  right-hand  side  of  (84)  can  be  written  in  a  somewhat  simpler  form  by  introducing 
the  function 


f  3  =  ■4i_ 

3  C*fn 


(85) 


and  observing  that  the  first  two  terns  within  the  brackets  can  be  written  as 


tof0(f3-l)  +^[f0(f3-l>]  =  7  fa  [f0(f,-D]  - 


<1<P 


df. 


d<p 


=  7fft 


«a 

dtp 


Use  has  been  made  of  the  equation 


—2  df0 

«  *•  *  *“'•  =  0  ■ 


(86) 


*  An  exception  Is  represented  by  the  case  of  axial  oscillations  (s^  =  0)  and 
4>  .  in  which  case  (84)  reduces  to  the  hypergeometrlc  equation2,3. 


? 


linear  with 
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which  f0  satisfies.  The  last  term  within  the  brackets  can  also  be  expressed  through 
f3  ,  observing  that  with  the  help  of  (86)  we  get  from  (83) 


\<f>  /  p<£  ' 

h f  (r)  ■  i\' -  '•  [  t 


±L.to^s\  =  £>(?,,  .mi 

icu  V  1  0  2  /  ia>  V  2  , 


Hence  (84)  can  be  written  in  the  form 


L($)  =  -  c2f0(ClF<1)  +  crp(2))  ,  (88) 

where  L($)  represents  the  left-hand  side  of  Equation  (84)  and  the  functions  F  are 
expressed  in  terms  of  f3  alone  by 


>(D  -  “vh 


— 


d^3  si/h  PQ  /_  ?  . 


pq  ;  F(2)  =  q2  — f+^T - (2czf  +  q£h  -  q2) 

d 4>  2La>  3  ta 


The  flow  rotation  can  also  be  calculated  from  (52).  Because  of  (55)  the  ^-component 
vanishes,  while  the  other  two  components  are  proportional  to  0  -  V1  .  In  view  of 
(79)  and  (82)  this  can  be  written  as 

..  ,  df„  dfn 

$'  -  V'  =  o*  — ~  +  C.  — -  . 

d0  1  d0 

From  the  third  expression  (87)  for  f2  and  (81)  one  obtains 

df 2  _  J_^£o  (-75 

d0  "  2i*>  d0  U  ‘ 

where  q^  =  [2/(7  +  l)]‘  is  the  value  at  the  throat,  obtained  from  (13)  by  setting 
q  =  c  .  Hence,  making  use  of  (86),  one  has 

-  V'  =  fcr(?  -  q^)  -  2i<yCj]  . 


and  the  vorticity*  is  given  by 


_ 2  1  v  L  v|  Vra 

*  Taking  into  account  that  the  variables  are  different  we  have 

V  *  dW  *  (8t,h/2[Wr]i)J„(sl,h[<’/<fiw]h  . 
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Observe  that  both  ?0'/r  and  r'i'Q  vanish  at  r  =  0  for  v  £  1  .  Only  for  v  =  1 
are  they  finite  at  r  =  0  and  only  for  this  v  can  there  exist  a  vorticity  on  the 
nozzle  axis.  Observe  also  that  the  subtractive  term  in  the  exponential  represents  the 
travel  tine  of  a  gas  particle  to  any  station  <p  ;  therefore  the  exponential  exhibits 
the  transport  of  vorticity  with  the  fluid  elements.  Clearly  for  isentropic  oscillations 
(a  =  0)  and  Cj  =  0  the  oscillations  are  also  irrotational,  but,  if  there  is  an  entropy 
perturbation,  vorticity  must  be  present  no  matter  how  Cj  is  chosen,  except  for  axial 
oscillation,  when  f'  -  0'  =  0  . 


9.  REDUCTION  OF  THE  SYSTEM:  BIDIMENSIONAL  NOZZLE 

The  developments  of  Section  8  can  be  repeated  almost  identically  for  the  bidimensional 
nozzle.  The  case  of  purely  axial  oscillations  (obtained  for  M  =  N  =  0)  provides  equa¬ 
tions  which  coincide  with  those  of  the  axisymmetric  nozzle.  The  difference  between 
the  two  only  appears  when  transverse  oscillations  are  present.  If  M  or  N  is  zero 
the  respective  equation  (64)  or  (65)  is  satisfied  identically  and  the  corresponding 
dependent  variable  Vv  or  Wt  disappears  from  the  equations.  All  these  cases  can  be 
obtained  as  particular  cases  from  the  case  when  neither  M  nor  N  vanish;  in  which 
case  we  obtain,  from  (64)  and  (65), 

Wj  -  V,  =  C0f0  ,  (90) 

C0  being  an  integration  constant  which  (contrary  to  the  case  of  circular  section) 
can  be  taken  non-vanishing. 

From  (66)  and  (67)  we  get 


a  =  s(0)  being  again  an  arbitrary 
Taking 


one  obtains,  from  (63), 

+  q7^'  +  P  =  o-f  t  .  (93) 

Combining  this  with  (64)  we  obtain,  as  before, 

-  V,  =  crf2  +  Clf0  •  (94) 

Hence  we  obtain,  from  (91)  and  (93), 


P,  -  czR, 

— * - *■  =  erf 

^  10  * 


(81) 


integration  constant. 


d$, 

U.  =  — 1 
1  i<t> 


(92) 


+  c^Rj  = 


i r  q2$J 


<r(fi  -Pf0) 
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and,  from  (74),  making  use  of  (90),  (92),  and  (94), 

toRj  +  Q2Rj  +  q5^'  -  (s^P2q2  +  cr£)$  =  “  (af2  +  cif0^s»°2c>2  +  +  CTNCof o  ■ 

Eliminating  R  from  the  last  two  equations  we  obtain 


q2(c2  -  ?)$"  -  7  (4  ~  +  2  to)  $! 

1  Vc5  d 4>  '  1 


2  7-1  toq5  dq* 


co  - 


7T-  W*  +0>2 


=  -a? 


i<y 


2  d<£ 

-  '•) +  ?  35  (? '  '•)  +  <s^?  *  ^)f‘ 

-  C^(s£V  +  or2)f0  +  C0c2o^f0  ,  (95) 


which  can  also  be  written 


Li^i)  =  ~  ^f0(C0P(10)  +  ciP,(1)  +(rp[2))  , 

where  Lj($j)  represents  the  left-hand  side  of  (95)  and 


(96) 


p(  0 )  —  ~2  . 

~  ~at,  • 


f!1^  =  s^p2q2  +  ; 


p‘!>  ■  «■$  ♦<*’<*♦■!> 


~T  "“y  ~T 

2N  2c  f3  +  qth  -  q 


2ico 


(97) 


From  (68)  we  obtain  the  vorticity.  Following  a  line  similar  to  that  of  Section  8 
we  obtain 


x  k  = 


c„b lT, 

v  2q 


Vevfr  - 


bp  —  ~  -  2toC- 

2  1  y 


exp 


ico  t  - 


d<£' 

V 


tb 


Again  we  notice  the  exponential  expressing  the  vorticity  property  of  being  transported 
with  the  fluid,  and  again  we  notice  that  the  vorticity  cannot  vanish  identically  with 
noniseutropic  oscillations,  except  for  purely  axial  oscillations,  when  =  Z'  =  0  . 


22 


10.  ADMITTANCE  CONDITION  AT  .HE  ENTRANCE  OP 
AN  AXISYMMETRIC  NOZZLE 


The  general  solution  of  (88)  is  made  of  the  general  solution  of  the  homogeneous 
equation  and  of  a  particular  solution  of  (88).  In  other  words,  if  we  know  the  solutions 
of  the  equations 


L(4h)  =  o 

(98) 

L($(j))  =  -c2f0P(J) 

(j  =  1.2)  . 

(99) 

the  general  solution  is 

$  =  Cj$<l)  +  o$<2>  +  C2$h  +  C3$h  .  (100) 

where  are  two  independent  solutions  of  (98)  and  C2  ,  C3  arbitrary  constants. 

Now  observe  that  (98)  has  the  singular  points  q  =  0  ,  q  =  c  =  q^  =  {2/(y  +  l)}»  and 
$  =  oo  .  For  a  supercritical  nozzle  with  a  finite  entrance  section,  only  the  sonic 
singularity  at  the  throat  is  of  concern.  If  the  solution  must  be  regular  at  the  throat, 
only  the  solution,  say  ,  which  is  regular  there  can  appear  in  (100);  and  the  parti¬ 
cular  solutions  must  also  be  regular  at  the  sonic  point.  Since  all  the  singu¬ 

larity  appears  now  in  the  independent  solution  ,  the  condition  of  regularity  at 
the  throat  is  simply  expressed  by  C3  =  0  . 

Hence  the  proper  solution  of  (88)  is 


$  =  Cj$(1)  ‘•o$<2)  +  C2$h  , 

with  ,  $(2)  ,  and  regular  at  the  sonic  point. 

From  (101)  one  obtains,  making  use  of  (79),  (82),  (55),  (80),  and  (78),  the 
expressions 

d$<n  d$<*>  d* 

u  =  c. — —  +cr — —  +  c,  — 

1  d  4>  d  0  2  d  <P 

V  =  W  =  Cj($'l,-f0)  +cr($(2>  -  f2)  +  C2$», 

P  +  q^U  =  -  C1iw$(n  -  cr(ia;§(2)  -  ft)  -  C2i«w$h 

S  =  Orf0  . 


(101) 


(102) 


If  this  is  considered  to  be  a  system  of  linear  equations  for  the  determination  of 
Cx  ,  cr  ,  and  C2  with  given  values  of  U  ,  V  =  W  ,  P  ,  and  S  for  a  given  <j>  ,  we 
see  that  only  three  of  these  values  can  be  arbitrarily  prescribed.  In  other  words,  a 
relation  must  ex. st  between  any  such  set  of  four  values,  a  relation  that  can  be  obtained 
from  the  relation  of  compatibility  of  the  four  linear  equations  (102),  that  is  from 
the  vanishing  of  the  determinant 


which  represents  two  aeiaittance  conditions  when  applied  to  the  entrance  of  tbs  nozzle. 
Here  &  represents  the  irrotational  admittance  coefficient. 

For  purely  axial  oscillations,  since  (48)  and  (49)  are  identically  satisfied,  one 
has  to  disregard  Equation  (82),  which  Bakes  use  of  (48).  The  corresponding  equations 
can  be  obtained  from  (102)  disregarding  the  second  one,  and  all  terns  in  C.  . 

c 
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With  the  same  procedure  used  above  one  finds  from  these  equations  the  corresponding 
admittance  condition 

U(<?£  +  M  +  PC  +  Sc5(io^:< 2>  -  f30  =  0  .  (106) 

For  isentropic  oscillations  this  condition  can  be  used  in  the  form 

£  a 

U  =  -t; - P  =  -  P  =  OCR  .  (107) 

q2£  +  ito  y 

In  this  case  a  corresponds  to  the  admittance  ratio  of  the  acousticians.  Observe 
that,  contrary  to  the  a  of  (i05),  this  a  has  to  be  calculated  for  st.h  =  0  . 

We  see  that  in  every  case  the  knowledge  of  three  functions,  £  ,  £(l)  ,  and  f(2) 
defined  in  (103)  is  sufficient  to  determine  the  admittance  condition.  The  numerical 
procedure  to  calculate  these  functions  is  explained  in  Section  14.  Here,  however,  it 
is  useful  to  obtain  a  more  explicit  expression  for  the  .  From  (98)  and  .99), 

recalling  that  L($)  is  given  by  the  first  member  of  (84),  we  obtain 

4  {*(J)L(V  -  *hL($(J>)}  =  (P  -  P)—r(P f0V(J)>  -  +  2H 

q^  d  <P  \c2  dp  y 


-2-$hF(J)  . 
>  n 


This  first  order  equs'.ion  has  the  general  solution 


?f0Vr"  =  e0l<t>)  constant  +  =v~=~  &==-  $hF(J)e~^(0>)  &4>'  ,  (109) 

q  (c2  -  q2) 


m  =  21"/  JTTj?  -l0‘e  (!-|?) 


and  4>'  is  again  an  integration  variable,  of  which  all  quantities  in  the  integrands 
are  functions.  The  lower  limit  in  the  integral  of  (110)  is  essential.  The  presence 
of  the  second  term  on  the  right-hand  side  of  (110)  causes  exp  0 3(P))  to  become 
Infinite  at  the  sonic  point.  Clearly  the  only  change  for  to  remain  finite 

there  (which  is  necessary  if  the  solution  must  be  regular)  is  the  vanishing  of  the 
integration  constant  in  (109),  since  the  lower  limit  of  the  corresponding  integral 
has  been  taken  at  the  throat.  Hence  the  proper  form  for  the  solution  of  (108)  is 


,(J>  = 


r/laj  2iw  \ 


bxp]\  [p  +  p-p 


F(J  >(<£') 

P(<P') 
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where  f0  has  been  replaced  by  its  expression  (77)  and  by 


immediately  obtained  from  (103). 

We  see  that  can.  in  principle,  be  obtained  from  (111)  once  £  is  known,  using 

only  integration  processes,  and  that  the  solution  of  the  Equations  (59)  is  not  necessary 
if  our  scope  is  limited  to  the  admittance  condition.  However,  in  practice  (111)  presents 
nasty  problems  of  evaluation,  since  -  as  the  form  of  the  exponents  shows  -  both  the 
integrand  and  the  factor  in  front  of  it  tend  tc  oscillate  faster  and  faster  when  the 
sonic  point  is  approached,  or  when  becomes  small,  toward  the  nozzle  entrance. 

These  oscillations  are  artificial  and  tend  to  cancel  each  other,  but  they  make  it 
difficult  to  carry  on  the  last  integration  which  would  require  absurdly  small  integra¬ 
tion  intervals.  A  better  way  of  calculating  the  *  is  to  make  use  of  the  equation 


—  [(?  - 
d<p 


[(c5  -  q5)^-*3] 


=  P(J)  ,  (112) 

a* 


the  validity  of  which  can  te  directly  checked  from  (111). 


11.  ADMITTANCE  CONDITION  AT  THE  ENTRANCE  OF  A 
BIDIMENSIONAL  NOZZLE 

For  the  bidimensional  nozzle  one  has  to  consider  one  constant  and  one  equation  more 
than  in  the  previous  case.  That  is,  if  the  solutions  of  the  equations 

Li($lh)  =  0  (113) 

Li($J3>)  =  (j  =  0,1.2)  (114) 

are  known,  the  general  solution  of  (96)  remaining  regular  at  q  =  c  =  (2 /(y  +  1)}'  =  q^j 
is 

=  C0$(I0)  +  +  o"£[2)  f  c2$,h  . 


where  all  partial  solutions  must  be  regular  at  qtb  . 

Hence  we  obtain,  from  (92),  (94),  (90),  (93),  and  (91), 

d$!2)  d$ 


lb 


d$(e)  d$‘n 

C„  — i —  +  c.  — - —  +  cr  — 1 —  +  C- - 

d<f>  1  d <p  2  d 4> 


v,  =  C0<k[0)  +  C, ($j 1  >  -  f0)  +  cr($<2>  -  fj)  +  C2$Ih 


vi-w, 


Vo 


(115) 
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+  0*0. 


=  -C0i<o$<0)  -  +  o-(fj  -  ico${2))  -  c2iw$ 


lh 


(115) 


S 


l 


<rtK 


Again  the  compatibility  condition  of  these  five  equations  between  the  four  constants, 
requires  the  vanishing  of  the  corresponding  determinant.  After  development,  and  divi¬ 
sion  by  f2$h  .  the  condition  becomes 


ujo1^”  -  0  -  i .«]  +  Vjtoc^*1*  -  i[0))  +  Wlto?#{°>  + 


+  p^J11  -  0  -  sx?  --(1)  * 


I1 '  +  -  f3cj  .  die) 


with 


1  d$ 


lk  . 


$ih 


£<J  >(<*>) 


c2f„* 


o  ih  L 


d 4>  lh  d<£ 


(117) 


When  applied  at  the  nozzle  entrance,  (116)  represents  the  admittance  condition. 

We  notice  that  if  the  axial  component  of  the  vorticity  vanishes  (Vli  =  Vt)  Equation 
(116)  becomes  identical  to  (104).  However,  the  functions  appearing  in  the  coefficients 
are  different,  depending  on  the  solution  of  different  equations.  The  isentropic  case 
again  can  be  derived  inuediately  from  (116)  taking  =  0  ,  and  the  isentropic, 
irrotational  case  leads  to  two  admittance  conditions  identical  to  (105),  the  only 
difference  being  the  subscript  1  on  the  various  quantities.  Finally,  the  case  of 
axial  oscillations  leads  to  equations  identical  with  (106)  and  '107).  Evidently  in 
this  case  there  is  no  difference  between  the  quantities  with  the  subscript  1  and  those 
without  subscript;  indeed  Equations  (84)  and  (95)  coincide  when  s.  .  =  s  =cr  =  0  . 
About  the  functions  and  ,  which  are  the  only  ones  to  be  determined  in  every 

case  for  the  purpose  of  obtaining  the  admittance  condition,  the  same  observations 
developed  at  the  end  of  the  preceding  section  for  the  corresponding  quantities  £  and 

1  apply  with  identical  conclusions. 


12.  SIMILARITY  OF  NOZZLES:  VELOCITY  DISTRIBUTION 
FOR  REFERENCE  NOZZLE 

The  solution  of  Equation  (88),  or  of  (98)  and  (99),  depends  on  the  parameters  u> 
and  svll  (</<w  has  a  fixed  value  once  the  reference  length  is  chosen;  see  footnote, 
p.  13);  it  depends  moreover  on  the  nozzle  geometry  through  the  function  q(<£)  (from 
which  c ('•/>)  and  p{4>)  are  derived).  The  results  obtained  for  a  given  nozzle  that  we 
shall  call  the  reference  nozzle  can,  however,  be  immediately  applied  to  a  whole 
family  of  nozzles  obtained  by  linear  deformation  of  the  axial  scale.  If  z  is  the 
axial  length  coordinate  measured  from  the  throat  and  positive  in  the  flow  direction, 
fi  is  tne  proper  scale  factor  and  qr(z)  is  the  axial  distribution  of  velocity  for 
the  reference  nozzle,  q(z)  =  (fiz)  will  be  the  general  velocity  distribution 
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for  the  derived  family  of  nozzles.  Prom  (22),  since,  in  view  of  the  one-dimensional 
assumption,  Ss  ?  dx  ,  we  obtain,  indicating  with  <£re{  the  potential  for  the  reference 
nozzle, 

0&<P  =  /?  q(z)dz  =  qref(/2z)dC#z)  =  d0ref  . 

It  appears  that  0{<t>  -  <£th)  =  4>ret  is  also  a  function  only  of  0z  =  z,.ef  ,  and  that 
as  a  result  q  can  be  considered  to  be  the  same  function  of  0ref  for  the  whole  family. 

It  is  immediately  seen  that  if  L (</>)  ,  that  is  the  first  member  of  (84),  is  divided 
by  02  ,  its  form  remains  unaltered  if  the  independent  variable  is  taken  to  be  <£_ef 
and  the  parameters  a >  and  s,,h  are  replaced  by  wref  =  on/0  and  s„h  ref  =  s vt/p  ; 
so  that 

(0 .  01 '  ref  ^ref  *  ^ref  •  svhreP  ’  (118) 

where  ^ref  *s  the  s°lutl°n  °f  (88)  obtained  for  the  reference  nozzle.  It  appears 
from  (103)  chat 


1  1  d$  1  d$ 

—  C(0 :  «*> .  svh)  -  =  ^ref^ref  -  "ref  •  s^h  ref>  ‘  (119) 

We  conclude  from  (118)  or  (119)  that  it  is  sufficient  to  solve  (98'  only  for  the 
reference  nozzle.  However,  the  solution  should  be  computed  not  only  for  the  discrete 
set  of  values  s„h  of  Table  I  (p.84),  but  for  s^  continuously  variable  in  a  certain 

range,  in  order  to  cover  the  possible  range  of  values  of  the  scale  factor  0  .  We  see 

also  that  if  the  nozzle  entrance  velocity  is  the  same  for  all  the  nozzles  of  the  given 
family  (that  is,  if  the  area  contraction  ratio  of  these  nozzles  is  prescribed),  the 
value  of  <pret  corresponding  to  the  entrance  section  is  also  the  same  so  that  £ref 
is  only  a  function  of  o>rej  and  svh  ref  .  For  given  scale  factor  0  and  mode  (svh)  . 

£ref  ,  as  well  as  £  ,  is  only  a  function  of  co  . 

Similar  observations  can  be  made  on  the  quantities  .  Inspection  of  (89) 

shows  that 

P(1)(0;sh)  =  0:Ff.*l(<fi  f  ;  svh  rpf)  ; 


F(2)  (<P\ o),  s„h)  =  0^Tl\(fPvet ;  "ref ;  s^h  ref) 


Hence  we  obtain  at  once  from  (111)  that 


1 5  ,  svtJ !  =  0f$(*nt  ;  *>ref  .  svl.  )  ; 

(120) 

§  (<t> ,  <o ,  s,,jj)  -  ^ret^ret  '  ^ref  •  svh  ref^  ■ 

We  are  particularly  interested  in  the  repercussion  of  these  relations  on  the  admittance 
condition  (’04) .  We  can  rewrite  this  condition  in  the  form 


U+AP+BV+CS  =  0  , 
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where  the  admittance  coefficients  are 

^(1)  -4 

A  "  -  O-  ico 

ta)pg(l) _ 

B  =  7(^<n  -  0  -ico 

■  +  !*,£<*>  -f 3}j 

C  =  ^(^<n  -  4)  - 


(122) 


(123) 


(124) 


Because  of  (119)  and  (120)  we  see  that 


A(0;  cu.  syh) 
B(0;  «.  svh) 
c(0 ;  o) ,  s„h) 


^ret^ret  ’  wref  ’  s^h  reP 
^®ref  ^ref  ’  ^ref  '  ref^ 
cref  ^ref  •  wref  •  svh  reP 


(125a) 


Hence,  from  the  admittance  coefficients  calculated  for  the  reference  nozzle  in  the 
appropriate  range  of  the  independent  variable  (or  of  q)  and  of  the  parameters,  one 
can  obtain  the  admittance  coefficients  for  any  nozzle  of  the  family. 


The  case  of  the  two-dimensional  nozzle  presents  the  same  property  of  nozzle 
similarity.  In  fact,  it  is  easily  seen  from  Equations  (113)  and  (117)  that,  defining 
s„  ref  =  s„//3  and  crN  ref  =  a^/P  ,  we  obtain 

^lh^  ’  w  •  sn  ’  ^n)  -  ^ih  ref  ^ref  ’  wref  ’  SK  ref  *  °"n  ref^ 

4^0;  8„  ,crN)  =  ref  (0ref  ;  0)ret ,  s„ ref  -  °"N  reP 

and.  from  (97)  and  (111), 

£ (i0)  (0 :  w  •  sM  .  <rH)  =  ^VefftVef  '<  "ref  .  s„  ref  .  ref) 
f(11)(0:^,  s„,orN)  =  ,^lVef(0ref:"ref  s*ref  •  CTN  ref> 

^[?)(0;c<;,  sM,crN)  =  ^IVef^ref  :wref  -  sgref  >0rNref)  • 


Here  again  the  subscript  ref  denotes  solutions  relative  to  the  reference  nozzle.  As 
a  consequence  if  the  admittance  condition  (116)  is  written  in  the  form 

U!  +  A,Pi  +B,Vt  +C1SJ  +  DjW,  =  0  . 


2b 


where  the  expression  for  admittance  coefficients  At  ,  Bt  .  Ct  ,  Dj 
obtained  by  comparison  witn  (116).  Then  it  follows  that 

A,(0;  w.  s„  .  o-N)  -  Aj  ref  (<^re/> .  wref  •  SM  ref  ,  crN  ref) 
Bj(0  ;  co  ,  sM  ,  o-N)  =  ^Bi  ref  ^ref  ■  wref  ’  SM  ref  ’  ref) 


is  immediately 
1 


(125b) 


C,(0;6),  s„.<rM)  =  Clref(^ref  ;  wref  ,  sMref  ,crNref) 

Dl(0;o>.  sM.orN)  =  /?D.  ref(0ref;  (orBt.  sMref  ,crNref)  . 

j 

Hence  from  the  admittance  coefficients  for  the  reference  nozzle  calculated  in  a 
sufficient  range  of  variation  of  the  independent  variable  and  of  the  parameters  one 
can  obtain  directly  those  for  any  nozzle  of  the  family. 

It  is  interesting  to  notice  that,  as  (125a)  and  ( 125b)  show  clearly,  the  scale  change 
from  the  reference  nozzle  to  any  nozzle  of  the  family  affects  in  a  different  fashion 
the  relation  between  scalar  or  axial  quantities  (pressure,  entropy,  and  axial  velocity 
components)  and  those  including  the  transveise  components  of  velocity. 


13.  NONLINEAR  ANALYSIS 

The  linearized  analysis  which  has  been  performed  applies  to  small-amplitude 
oscillations  and  is  most  useful  in  the  treatment  of  spontaneous  inst ibilities.  It  can 
be  used  in  the  prediction  of  the  stability  of  the  steaiiy-state  nozzle  operation  and, 
if  the  regime  oscillation  in  the  unstable  situation  has  a  small  amplitude,  it  can  be 
used  to  predict  some  characteristics  of  the  oscillation.  However,  if  the  oscillation 
does  not  initiate  spontaneously  but  instead  requires  a  finite-size  disturbance  to  the 
steady-state  operation  in  order  to  excite  an  oscillation,  the  linearized  analysis  is 
not  sufficient.  Also,  if  the  regime  oscillation  does  not  have  a  small  amplitude  (as 
is  often  the  case),  the  linearized  analysis  does  not  accurately  predict  all  of  the 
characteristics  of  the  oscillations.  In  these  situations,  a  nonlinear  analysis  is 
better  suited  on  the  basis  of  accuracy. 

The  analysis  of  the  axisymmetric  nozzle  was  extended  to  include  nonlinear  effects 
by  Zinn8.  A  perturbation  series  was  employed  where  the  perturbation  quantity  was  an 
amplitude  parameter.  Of  course,  the  first  order  solution  is  identical  to  the  linearized 
solution  discussed  in  previous  sections.  The  second  and  higher  order  solutions  re¬ 
present  the  nonlinear  effects.  In  his  work,  Zinn  completed  the  analysis  up  to  and 
including  third  order;  however,  for  the  sake  of  brevity  only  the  second  order  results 
will  be  presented  here  and  the  reader  is  referred  to  the  original  work  for  third  order 
results. 

The  same  nondimensional  scheme  given  in  Section  2  has  been  employed  so  that  the 
nonlinear  system  of  equations  describing  the  oscillatory  phenomenon  is  given  by 
Equations  (5)  through  (9).  It  is  convenient  to  transform  the  time  variable  by  setting 
y  =cot  ,  where  cu  is  the  angular  frequency,  so  that  the  solution  is  277-periodic  in 
y  .  In  accordance  with  the  standard  procedure  in  nonlinear  mechanics,  the  following 
series  are  substituted  into  the  nonlinear  system  of  equations: 


q  +  £q{1)  +  £2q(2)  +  0(£3) 
p  +  ep(  15  +  e2p<2)  +  o(£3) 


p  -  p  +  ep(1>  +  e2p(2)  +  o(e3) 
s  =  8+es(1)+e28(2)+  0(e3) 
go  =-  coi0)  +  ew(1)  +  e2w(2>  +  0(e3)  , 

where  £  is  the  amplitude  parameter  and  barred  quantities,  as  before,  denote  steady- 
state  solutions.  Separation  of  the  equations  according  to  powers  of  £  yields  the 
same  equations  for  the  steady-state  quantities  as  (11),  (12),  and  (13).  The  equations 
for  the  first  order  coefficients  q(1)  ,  p(1)  ,  p(1)  ,  and  s(1)  are  the  same  as 
Equations  (17)  through  (21)  for  q'  .  p'  ,  p'  ,  and  s'  ,  except  that  the  partial- 
differential  operator  d(  )/3t  is  replaced  by  co(0)B(  )/3y  . 

The  second  order  system  of  equations  becomes 

Continuity: 


3p<2> 


+  V.(qp<2)  +  q(2>p)  =  -  V«  (p(1>q(1>)  -  oo(l) 


3p(1) 


Momentum: 


3q<2) 


+  ^7(q.q<2>)  +  (V  x  q(2))  x  q  + - V(q.q)  +— rVp(2) 


f0,  P(l)3q(1>  p{1) 

=  -  or  ' -  +  co- u '  — - +  — =— 


+  —z —  V  (q.q*  *>)  + 
p  By  p 


1 

+  1  —  -  (V  x  q(l))  x  q  +  -  V(q(1>.q(1>)  t  (Vx  q111)  x  q(1) 
P  2 


Entropy  Equation: 


3s<2)  3s(1) 

a><°>  —  +  q.Vs(2)  =  -  q^.Vs*1’  -oo<l)^— 
By  ay 


Equation  of  State: 


(2)  P(2'  P(2) 

r  ’ - 3-  +  _ 

7P  P 


1  If  Ell 

2  y  \  p 
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As  done  In  the  linear  analysis,  a  scream  function  and  a  potential  are  defined  for 
the  axisymmetric  case  according  to  Equations  (14)  and  (16).  Furthermore,  they  are 
still  used  as  independent  variables,  as  discussed  in  Section  4.  Then  we  have  (con¬ 
sistent  with  (27))  u  ,  v  ,  and  w  such  that  q  =  ue^  +  ve^  +  wee  .  Now  the  equations 
for  p(1)  ,  p(1)  ,  s(l>  ,  u(1)  ,  v(1)  ,  and  w(1)  become  identical  to  Equations  (28) 
through  (33)  for  the  primed  quantities  except  that  again  the  time  derivative  has  been 
transformed.  The  equations  for  the  second  order  quantities  become  the  following: 

Continuity: 
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'//-Component  of  Momentum: 


v<2)\  [  ^  ^  /v(2\  ^  fv(2\  _  ^(l)  ^  / v<1^j 

pqrj  3</>\PQr/  3^  KtP J  3y\pqr/ 


F <r  %  f2/ v(1) 

2  3</.‘  r  Vpqr 


1  3  (rw(1))  u(l)  B  A(l  j  p{1)  ^  (P<1  )  - 

2pqB^  r  ^  q  3</>\pqr/  p  B '/'\?/P/  B'/' 


^-Component  of  Momentum: 


r)  /  (2y 

w(°)J_(rw(2))  +^A(rw(*))  4  ?^ 


B0  \yp 


B  B/v(1)V 

By  30\yjqry 


#39  «  3<*>  p  39\7P/  39 


Entropy: 


■B8<2)  _  as<2> 

„co)  « —  + 

By  3</> 


=  -co<l> 


-gu(1)B8(n 
Q  q  B^ 


Equation  of  State: 


,  v<l)  3s(1)  pq  fn  3s<l) 

-  2p#ii;~wr  3— 


TO)  If"'  ■>'”  i 

p  yp  2  v  p 


y  A(1)N\ 

c7  \7P/ 


=  F  . 


(134) 


(135) 


The  first  order  condition  given  by  (52)  and  (55)  (zero  axial  vorticity)  has  been  used 
in  deriving  the  above  equations. 

The  time  dependence  of  the  first  order  solutions  as  shown  by* (44)  and  the  form  of 
the  inhomogeneous  parts  of  (130)  through  (135)  indicate  that  the  y  dependence  of 
those  inhomogeneous  terms  is  readily  expressed  as  a  Fourier  series.  This  indicates 
that  second  order  solutions  are  of  the  form 

u(2) 

-  =  5Z  elny£  (<£,</»,  0) 

q  n 

v<2) 

—  =  51  eInV0.  ./,.£?) 

pqr  n  ^  i  (136) 


«'  -  -ir-JF-t'J 


•  ^  ***** 


rw<2)  =  2H  einy£n(<£,i p ,  6)  [ 


->(  2 ) 


=  £  elny*n(0,  0,  (9) 


P  n 

,<2> 


—  =  2Z  einywn(0, 0 .  0) 


7P 
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(136) 


j(z)  -  5Z  ela*o-  (<fi'\p,&)  . 


As  shown  later,  the  solutions  contain  only  terns  corresponding  to  n  =  0  and  n  =  2  . 

Substitution  of  (136)  into  (130)  through  (135)  and  separation  of  the  Fourier  com¬ 
ponents  yields 


W°>k„+?  i^  +  ^^  +  2pqjL(t/,^)  +^ia  =  ^ 


30  30 


30 


2  0  20 


,  (Q)C  ^  ,-y-  v  1  dQ7  3b. 

in<yl  J,0  + —  (<r£_)  +  — a - cr  -  =  — a 

n  30  n  30  2  “  d0  30 


lncn' 


J0)r 


^n  _ 

** 

30  30 

30 

-,3L 

a  — **■  +  — = 

lEa 

30  30 

30 

“■w  =  ^ 

C*  K  +  ~  Pn  - 


(130a) 


(131a) 


(132a) 


(133a) 


(134a) 


(135a) 


where  ^  .  31^/30  ,  etc. ,  are  the  coefficients  in  the  Fourier  expansion  of  A  , 
3B^3£  ,  etc. 


In  analogy  with  (7S)  we  define 


<f  s 
n  30 


(137) 


Then  (131a),  (132a),  (133a),  and  (137)  yield,  upon  combination  and  integration, 

t<t> 

I 

'o*f. 


^  =  f°n30 


E&_LE&d0'  +  %.  -0n(c itb,0.P) 

*  Aw» 


th 


k0Q 


"on 


+  fonV^b  ’*•  &  (138> 
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=  fon35 


L’0th 


Pg_...  .G°  44>>  +  i  _  $  (0  ,0,0) 
«  fon  fon 


+  fo7n^th  •  •  <139> 


where  the  definitions  are  made  that 


if*  d7 

Uth 


*on  "  e 


•  in&) 


«»  r  d0' 


7  ' 


th 


In  addition,  (134a)  is  readily  integrated  to  obtain 

r '4> 


CT  f 

n  Acn 


i*. 


•=^-d<^'  +<rn(<t>th  0) 


th 


■■on 


(140) 


Using  (131a),  (135a),  and  (137)  through  (140)  to  substitute  into  (130a),  we  find  that 


3^ 


+  + 


-  7  2W°>  .  . 

B<£2  y  c*  d <ty 


nw 


(0) 


?  2 


<#> 


+ 


+  2Pq7—  (</> 


3$, 


Bi/;  \  2</> 


s i  .fE  f*k  =  ?  C  t .  .  y-L/KisS  . 

2 </>  B52  2,1  \  e2  /] 


inw(0)(Fn  +  Gn)  -  2PQf, 


CD 


±L3J^ 

i  + 1  Vh» 

d'P  \  B 4>) 

1  A'P  302 

In  •  <141) 


where  Hg  is  defined  as 

"•f 

J*tb 


G°  d0'  +  ^?n(0th  .  'P .  6)  -  $D(<fitb  .  t .  5) 
"  ron 


and  7^  is  a  potential  function  defined  by  the  relation 


The  inhoaogeneous  part  IQ  in  (141)  has  been  simplified  by  assuming  zero  axial  vorticity 
to  second  order. 
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The  form  of  Equations  (137)  through  (141)  and  the  form  of  the  first  order  solution 
(as  given  by  (44)),  together  with  the  boundary  condition  that  the  radial  velocity 
vanishes  at  the  nozzle  wall,  indicate  how  the  radial  and  tangential  dependencies  of 
the  second  order  solution  may  be  obtained;  the  solution  may  conveniently  be  found  in 
the  form  of  an  expansion  in  eigenfunctions.  Following  standard  procedures,  Zinn  obtains 
the  following  for  a  wave  traveling*  in  the  negative  d  direction 


while  for  a  standing  wave  he  obtains 


If  n/0,2.  I<s)  = 


=0  .  The  symbol 


Z 

i=  n  : 


means  that 


cos  w6 


m  =  0  and  m  = 


(142b) 

2 


produce  the  only  nonzero  values  of  In  mv  „  .  Note  that  <y(1)  was  properly  set  equal 
to  zero  in  the  above  relations  for  the' combustion  instability  problem  considered  in 


Reference  P.  Otherwise,  the  first  harmonic  would  appear  above.  Note  that  s  „ 
is  the  q11*  root  of  J„v(x)  =  0  .  Solutions  are  found  in  the  following  form  for 
traveling  waves 


Since  transverse  standing  waves  are  the  sua  of  two  transverse  traveling  waves  in  the  linear 
case,  no  distinction  was  necessary  there.  However,  here  in  the  nonlinear  case,  it  Is 
necessary  since  no  simple  relation  applies  between  the  twe  types  of  saves. 
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^2n  ”  *on 


f“£  -5 

1  dqz 

1  f  ll.  id)'  . 

2  on  <tf>'  V 


Integration  of  (131a)  and  combination  with  (137),  (143a),  and  (150)  determines 
Pn  UP  .  Then  this  result  is  combined  with  (135),  (143),  and  (150)  to  determine 
R„  ‘ ’  .  The  results  are 

□  ,  D*7  ,  Q 


1  dq 


n-m,/'Q  2  dtf>‘ 


7Sn.■v>Qd^>,  T  Bn,av.q  inw<  ^n.av, q  "  q  .^^n.av.q  ^151^ 


Rn.«l'.q  “  T5  ^n.mv.q  +  Pn,«K,q^  ~  Sn,mv.q 


Obviously,  since  the  coefficients  ,  Bn  %v  q  ,  etc.,  are  different  for 

traveling  and  standing,  so  will  the  solutions  given  by  (146)  through  (152)  be  different. 
Note  further  that  the  oj  in  the  linear  solutions  given  by  (102)  actually  represents 
the  w<0'  of  (144)  through  (152).  In  the  linearized  problem,  it  is  not  necessary 
to  distinguish  between  the  two,  since  the  difference  is  of  higher  order.  Similarity 
is  noticed  between  (104)  and  (147)  through  (151). 

Substitution  of  (146)  with  C3n  mv  Q  =  0  into  (147),  (148),  (150),  and  (151)  leads 
to  four  equations  which  are  linear’ in’the  th.ee  constants  crn  q  ,  C1D  By  Q  ,  and 
C2n  mv  Q  .  In  similar  fashion  to  the  development  in  Section  10, 'the  requirement  of 
compatibility  leads  to  the  following  relationship  between  Un  Bl/  ,  Pn  -l/  q  ,  Vn  uv  Q  , 

^n.av.q  • 


^n.av.q  • 

Un.mv.Q£?^n!av.q  ~  Sn,«v,q)-  inw(0)]  + 


+  V  i nu)^  3  ^  C  ^  ^  +  P  (C2£(n  —  £  )  — 

vn,av,qiD  c  sn,av,q  +  rn.av,q'c  sn,av,q  '’n,av,q/ 


q  ~  ^3n^n,av,  q 


=  W0)Pfon^fiVi(j  +  + 


~5  _  T 

+  v<2)  (£V(l>  -£  )  -  id3’  c5  q  ~  °^h  A<  n  + 

^n.av.q''1'  =n,av,q  ^n.av.q'  ''n.av.q''  2  ^n.mv.q 


+  inco(  3  ^  ^ 2  *  —  f  C 

“  ^n.av.q  x3nvi 


I . av , q 


*  *  ft „t  •  I-t*?5? *&° *  * 
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The  following  definitions  have  been  made  in  order  to  obtain  the  simplified  version 
(153) 

^n.mv.q  =  $  'JT^hn.mv.q  ^154* 

ha.rnv.q  us" 

f( J)  r—  _ i _ <£(J)  _£_  $  _  &  4  >f>(  J  )  I  /]CB\ 

bn.mv,q  “  jj5fo^^  IJn.mv.q  ^  xhn,mv,q  rhn,mv,q  ^  xn.mv,qJ  > 

where  J  =  1,2,3  ,  ^ 


v(  l)  =  f 

"n.mv.q  Aon 


.  _  p  —  R 

T?f  n.it'.q  °n,Bv#a 

**  1On 


f<t>'  /P*" 

2  is!  f  ( 

2  d <p"  on\ 

J*th  V*tb 


%*1U1 Ld0»/  W"  d<y 
^fon  / 


tf>  — j  /  t*'  \ 

(2)  =2  —  f  I  Sje,<Lfrh"  ]d0'  +  B 

‘■■"•a  2  d^  onl 


Ki3)  =  f„  -En--l^.:iL  d0'  . 
"n.wv.q  •‘•on  -ryf 


Equation  (153)  applies  everywhere  along  the  convergent  portion  of  the  nozzle  and, 
when  applied  at  the  nozzle  entrance,  is  the  nozzle  admittance  expression.  The  simi¬ 
larity  with  (104)  is  immediately  seen. 

In  the  irrotational  case  where  and  Cln  B  are  zero  the  nozzle  admittance 
expression  is  simplified.  In  that  case  (147),  (148),’  ana  (151)  yield,  after  application 
of  the  definitions  (154)  and  (155), 


0  _  C  V  =-c^f<?*3* 

wn,*i',q  =n,av,  q'n.av.q  J-on=n,av.q 


+  ina>(0)U, 


n.ai'.q 


+  ^n,av,q^n,av,q  +  0  ^n,mi\q^ 


T  irTzj  ina>(0,c2f  <f(3) 
:>n,a»',qRn,av,q  lnw  c  Ion=n,*»',q 


The  problem  remains  to  determine  the  functions  £n  tv  q  and  q  to  substitute 

into  the  relations  (153)  and  (156).  Reference  8  shows  that  Cn  1  is  governed  by 
a  complex  Riccati  equation  in  a  similar  manner  to  the  function  £  in  the  linear  case 


f 
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of  Section  10.  Also,  if  oj  were  replaced  by  n ro(0>  ,  £  by  bjr 

,  and  J)  h~  :  q  in  (112),  the  governing  equation  for  Za.mv.n  is 

oh'aiied.  So  *'■  :■  saia-i  ...  1  ;al  techniques  may  be  used  to  deteraine  CnBV>q 
gtj)  as  are  used  to  deteraine  £  and  .  These  aethods  are  discussed  wit 


^n^ay  as  used  Co  deteraine  £  and  .  These  aethods  are  discussed  with 

respect  to  the  linear  theory  in  the  next  section. 


There  is  one  iaportant  difference  between  the  linear  case  and  the  nonlinc.^'  case. 
This  difference  appears  on  account  of  that  inhomogeneous  part  of  (144)  proportional 
to  F^3^„  q  and,  consequently,  on  account  of  the  solution  of  the  function  ^wv.q  • 
F^33v  q  is  a  sum  of  terms  each  of  which  is  proportional  to  the  squares  or  double 
products  of  the  first  order  solutions  for  the  flow  properties  as  given  in  (102). 
Therefore,  in  order  to  calculate  the  second  admittance  coefficients,  it  is  necessary 


to  calculate  the  first  order  flow  properties. 


in  (102) 


should  be  known  to  completely  determine  the  second  order  admittance  coefficients.  That 
would  require  a  coupled  solution  of  the  oscillatory  phenomena  in  the  nozzle  and  in  the 
chamber  preceding  the  nozzle. 


There  is  a  prefeired  method  of  solution  which  uncouples  the  numerical  integration 
of  the  nozzle  equations  from  the  analysis  of  the  chamber  flow,  although  it  was  not 
employed  in  Reference  8.  Due  to  the  nature  of  F*3B„  q  and  of  (102),  F*3B„  q  and 

^n3*u/  q  nay  be  tbe  convenient  forms 


n.av.q 


S^)2  +  s2(c2)2  +  a3(or)2  +  S,(ClC2)  +  Ss(Cl0r)  +  S6(C2o-) 


^njwv.q  =  ai(Ci)2  +  a2(C2^2  +  a3^Z  +  a«(ClC2>  +  a^CjCr)  +  a6(C 2cr)  . 


are  independent  of  C. 


Since  the  coefficients  of 


R11C1 C  OUU  ai  C  AUUCyCU’.'CUU  ui  v.  ,  Vj  •  Miu  V  .  WlUVC  WUt  wv&aAAV/avuvs  V* 

2,  are  constant  throughout  the  nozzle,  may  replace  F*3gV  q  and  cci  may  replace 
^aav  i  in  tbe  ®odlfied  version  of  (112).  Of  course,  the  CCj*  s  may  be  determined 
by’ numerical  integration  without  a  knowledge  of  the  constants  Ct  ,  C2  ,  and  cr  . 

Note  that  Q  ,  K^2^  ^  ,  and  may  also  be  written  in  identical  form  to 

(157)  so  that’ the  inhomogeneous  parts  of  (153)  and  (156)  may  be  presented  in  that 
same  form,  where  Cj  ,  C2  .  and  cr  remain  to  be  determined. 


Note  that  there  are  actually  a  double  infinity  (two  n  values  and  an  infinity  of  q 
values)  of  relations  (153)  and  (156)  for  traveling  waves.  In  tbe  standing  wave  case, 
there  are  a  quadruple  infinity  (two  n  values,  two  m  values,  and  an  infinity  of  q 
values)  of  those  relations.  Each  of  these  infinite  number  of  relations  is  to  be  applied 
as  the  admittance  relation  for  a  particular  eigenfunction.  In  a  practical  sense,  con¬ 
vergence  will  occur  with  a  finite  number  of  terms. 

(fciite  similar  relations  to  (153)  and  (157)  are  presented  in  Reference  8  for  the 
third  order  (in  c)  quantities.  Of  course,  the  higher  tbe  order,  the  more  cumbersome 
are  the  algebraic  developments  and  the  numerical  integrations  to  be  performed. 


PART  II,  APPLICATIONS 


14.  CALCULATIONS  FOR  “CONICAL”  NOZZLES 

The  vast  Majority  of  nozzles  in  practical  devices  are  axisyanetric,  with  the 
generatrix  of  the  convergent  portion  (Pig.  17)  consisting  of  a  straight  segaen^ .  generating 
a  convergent  cone,  and  circular  arcs,  generating  the  throat  portion  and  possibly  the 
entrance  portion  (which  connects  the  conical  convergent  to  the  cylindrical  coabustion 
chaaber). 

The  calculations  were  perforaed  for  a  reference  nozzle  of  this  type,  using  the 
equations  developed  in  Part  I  for  the  axisyaaetric  case.  In  order  to  include  in  a 
single  calculation  nozzles  of  various  contraction  ratios,  the  entrance  portion  of  the 
nozzle  was  assuaed  to  be  part  of  the  conical  convergent,  so  that  the  slope  is  dis¬ 
continuous  at  the  entrance  (Fig. 17).  Within  the  one-diaensional  approxiaation,  this 
allows  choosing  each  station  in  the  calculations  to  be  the  nozzle  entrance  station. 

The  effect  of  the  presence  of  an  entrance  portion  other  than  conical  was  studied 
separately  (Section  15)  and  found  indeed  to  be  generally  saall. 

It  aust  be  observed  that,  according  to  Section  12,  a  whole  faaily  of  nozzles  can 
be  obtained  by  axially  stretching  or  shrinking  the  reference  nozzle.  However,  the 
generatrix  of  the  nozzles  thus  obtained  is  Bade  of  straight  segaents  and  elliptical 
(rather  than  circular)  arcs.  This  Beans  that  the  correspondence  between  a  practically 
designed  nozzle  (different  froa  the  reference  nozzle)  and  a  nozzle  of  the  faaily  is 
only  approxiaate  in  the  throat  (and  aaybe  the  entrance)  portions. 

In  Part  I  the  physical  quantities  have  been  nondisensionalized  with  respect  to  the 
stead7-state  stagnation  properties.  Indeed,  this  is  the  aost  appr>  priate  choice, 
because  the  saae  stagnation  properties  can  be  used  as  reference  quantities  in  the 
study  of  the  coabustion  chaaber  behavior,  which  Bakes  use  of  the  sdaittance  condition 
at  the  nozzle  entrance.  For  the  saae  reason  it  is  logical  to  choose  the  coabustion 
chaaber  radius  (or,  in  the  two-diaensional  case,  its  height)  to  be  the  appropriate 
reference  length  L*  introduced  in  Section  2. 


On  the  other  hand  then  this  study  was  undertaken  several  years  ago  the  reference 
quantities  were  taken  to  be  the  static  properties  at  the  throat  and  the  throat  radius 
(so  that  the  nondimensional  steady-state  pressure,  density,  velocity  and  radius  at 
the  throat  took  the  value  of  unity),  and  the  initial  programing  for  the  computer  was 
developed  on  this  basis.  In  subsequent  calculations,  rather  than  changing  the  program, 
it  was  decided  to  stick  with  this  choice  so  far  as  the  calculations  were  concerned, 
only  switching  to  the  stagnation  reference  properties  when  the  final  admittance  co¬ 
efficients  were  calculated. 

Although  this  may  be  a  little  confusing,  the  decision  was  made  to  keep  in  this 
Part  II  the  alternate  choice  of  the  reference  system  so  as  to  avoid  the  laborious 
conversion  of  numerical  results  or  analytical  developments  based  on  the  throat  refe¬ 
rence  system  to  the  stagnation  reference  system.  However,  we  shall  use  for  the 
variables  the  same  notations  used  in  Part  I,  explicitly  stating  when  the  stagnation 
reference  system  is  used  instead  of  the  throat  reference  system. 

In  view  of  the  relation 


'-LYr 

7*  V 


which  immediately  follows  from  (13)  and  from  the  definition  (9),  Equations  (12)  and 
(13)  become 


P  _  -7  -y.j  7  +  1  7-1-s 

—  -  cz  =  p'  - - q2 

P  2  2 


Of  course,  since  at  the  throat  we  have  p  =  p=  q=  c  =  rw  =  l  we  obtain,  from  (43) 
Equation  (88)  becomes 

L($)  =  q7 (q^  -  c5)#"  -  q2  (  2iw  +  — )  $'  + 

V  2c5  d  4>J 


+ r  ■  ia>  I  ■  *c2r/(r'l)  ■&)  * 


-  (CjF(:)  +  aF(2))PfQ 


where,  instead  of  (89),  (77)  and  (85),  we  have 


«Z„Q52/(T-,> 


Svh^C 


=  2/(7-:) 


(2c*f3  +  1  -  q5) | 
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The  appropriate  solution  of  the  equation  for  $  is  given  by  the  equation  corres¬ 
ponding  to  (101). 

$  =  Cj$(1)  +cr$(2)  r  C2$h  . 

where  all  the  &  s  Bust  be  regular  at  the  throat.  The  adnittance  relation  can  be 
written  in  the  fora  corresponding  to  (121), 

U  +  AP  +  BV  +  CS  =  0  .  (159) 

with  U  .  P  .  V  and  S  still  defined  by  (44)  and  A.B.C  bv  (122).  (123)  and  (124) 
with  q^  =  1  . 

This,  however,  is  not  the  aost  convenient  fora  for  use  as  a  boundary  condition  on 
the  oscillatory  flow  in  the  coabustion  chaaber.  Indeed,  the  equations  in  the  cylindri 
cal  chaaber  are  usually  separated  by  setting 

u'  =  U(z)J„(svhr)0(^)ei£,'t 

v'  =  U(z)  ^  (svhr)e(0)elat 
dr 

(160) 

p'  =  P(z)JJ/(s„>1r)6(^)ri1"t 
s'  =  i(z)Jv(svhr )melut  , 

where  r  is  nondijensionalized  by  the  chaaber  radius,  u'  ,  v'  by  the  stagnation 
sound  velocity  c°  and  p'  by  the  stagnation  pressure. 

Coaparing  (160)  with  (44),  and  taking  into  account  the  difference  of  the  reference 
systeas,  we  obtain 


After  substitution  the  adaittance  condition  becoaes 


U  +  dP  +  SU+Ci*  =  0  , 


(161) 
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where 

\(y  ♦  i\*<yo/(yn  q 

2  /  g2/(r-i)  q?(c^(I)  -  £)  -  icj 

* 

qi)  7%(1> 

B  "  1  ci/(y'l)  F(^(1>  -  0  -  iw 

-  /  2  V— 5  {(1  -  ?)£{1>  -i^?)  +  f,C 

\y  +  1/  ac  TO^** -o-iw 


Certain  (otherwise  redundant)  additional  definitions  are  aade  necessary  by  the  way 
the  original  coaputer  prograa  was  written.  Instead  of  <t>  we  introduce  a  new  indepen¬ 
dent  variable 

4>  =  2k{4>  -4>th)  ■ 

The  value  k  derives  froa  the  following  slaple  considerations.  Since,  by  definition, 
d <t>  -  5dz  .  we  find  d<£  =  2<qdz  .  with  z  representing  again  the  physical  axial  co¬ 
ordinate.  Hence 

dq^  __  1  dq 
x  dz 

The  particular  choice 

K 

results  in  dq^/d#  =1  at  =  0  .  This  is  the  value  of  k  that  was  adopted  .n  the 
following  calculations.  For  the  throat  geoaetry  of  Pigure  17.  x  can  be  expressed 
explicitly  (using,  for  instance,  the  procedure  outlined  in  Reference  9)  as 


In  conjunction  with  the  independent  variable  the  additional  definitions  were 

aade 

£  =  2xl  .  #(l>  =  2 x£(l)  ,  co  =  xco  ,  =  *svh  . 

In  teras  of  the  variables  the  adaittance  coefficients  aay  be  written  as  follows: 

i(y  +  Mcyt)/(y-t)  q 

U  =  y\  2  )  c2/(r':>  < 

Q  _  _ C?|(U _ 

lx5i/(y-i)  ?(g5|d>  _q  _ii(j 


(162) 


i»H,  Jf** *  *"  "f*1' 


2  v_. +  r,£ ' 


7  +  V 


q^c5!* n  -  6  -  i 


The  adaittance  coefficients  are  coaplex  nuabers  since  £  .  i(1>  .  gl2)  .  f3  and  i&> 
are  coaplex. 

An  interesting  observation  is  that  the  second  order  differential  equation  (158) 
need  not  be  solved  for  .  $tl)  and  $(2>  .  ^It  is  sufficient  indeed  to  solve 
directly  the  equations  which  are  satisfied  by  £  .  and  £(2>  .  that  is* 


(g  +  IbK  -0  -  i*> 


[(1  -  ?)|(l>]  -= 


r  ~  (—  2  \ 

\2<?  (7  +  i)(i  -  a7)) 

i 


-  w-j  * 


2 

+ - -2* -  (164) 

7  +  1  4  q 


-•  £  -  (~=f  +  - - - — •  [(i  -  q 5)f(2)J  + 

\2q2  (7  *  i)d  -  q2y 

7+1  &<P  2iw  V  2q2  q2  * 


where  the  following  definitions  apply: 


b  =  q^Cc7  -  q5) 


7  +  l?d? 


c?  d<£ 


h  =  coq5 


«2  qc^/O'-D 
1  — ; — *•* 

y  -  i  ?  a?  . 

k  =  - — — s-<a  . 

4  c7  d</> 

Equation  (163) ,1s  the  Rlccatl  equation  obtained  directly  froa  (158)  upoo  substitution  of  ( 
for  (l/{>)(u'?/dd').  Equations  (164)  and  (165)  are  obtained,  after  the  proper  adjust  went,  froa 
Equations  (112). 
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A 

The  nonlinear,  first  order  equation  for  £  is  a  complex  Riccati  equation  and  may  only 
be  solved  by  numerical  integration.  Once  this  is  done,  the  linear,  first  order  equations 
for  ari  £<2)  may  be  solved,  obtaining  the  standard  integral  forms  which  are 

essentially  identical  to  the  form  given  by  (111).  However,  rather  than  numerically 
evaluating  the  integral  solutions,  it  was  more  convenient  to  solve  all  three  complex 
(or  six  real)  equations  simultaneously  by  numerical  integration. 

The  numerical  integration  requires  statements  on  the  initial  conditions.  These 
are  provided  by  the  condition  of  regularity  at  the  throat  (4>  -  0).  After  assuming  a 
Taylor  series  expansion  about  the  throat  and  evaluating  the  coefficients  in  the  series, 
we  find 

£(0)  =  cc0  +  t(30 

=  0 

£(2)<0)  =  0  . 

where 

7  +  1  <S2  -  s2h  y  -  1  A2 
- - - or 


2  4  4 


Since  the  Riccati  equation  is  singular  at  the  throat  (c7  -  q*  =  0),  the  numerical 
integration  cannot  begin  exactly  there.  So,  in  addition  to  the  above  conditions  the 
first  derivative  of  t  at  the  initial  point  must  be  obtained  from  the  regular  expansion 
and  supplied  to  the  numerical  integration  scheme.  This  initial  condition  is 


where 


ctj  +  . 


a  s  _  (y  -t  l) 2 

1  co2  +  (y  +  1) 2  |  2 


(a2  -  0i)  +  (y  +  1  )£aJ0  + 


A  (7  + 
or  +  - 

4 


l)3 


+  (y  +  i)2S 


.]  .  (y  +  i)  ( 

- 4 
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fix  = 


1 


co2  +  (7  +  iy 


(7  +  D2a/0  -^|-^(a2  -  $)  + 


1-7  -K  ~ 

+  (7  +  1) - wa0  -  (7  +  i)wa0b  + 

2 


a?  + 


(7  +  D: 


7  -  1  „ 


co  [s2h  +  (7  +  I)2)  + 


72  -  1 

- wb 


Note  that^  b  is  a  coefficient  in  the  power  series  expansion  for  q7  .  That  is 
q2  =  1  +  4>  +  b<£2  +  ...  ,  the  value  of  b  being  determined  by  the  nozzle  geometry  in 
the  throat  vicinity.  The  above-mentioned  power  series  shows  the  convenience  of  the 
transformation  of  the  independent  variable  from  4>  to  <jt>  ;  the  coefficient  of  the 
first  order  term  in  the  power  series  is  unity. 


Rather  than  numerically  evaluating  the  integral  for  f3  ,  it  is  convenient  to  solve 
a  first  order  differential  equation  for  f3  simultaneously  with  the  equations  for 
£  ,  ,  and  £(2)  .  As  directly  checked,  the  proper  equation  is  as  follows: 


d 

d<£ 


Id? 
2  "d0~ 


(166) 


subject  to  the  initial  condition  that  f3(0)  =  0  . 

The  steady-state  velocity  profile  q{<£)  must  be  determined  for  tnc  given  geometry 
of  the  convergent  portion  of  the  nozzle.  This  cannot  be  found  in  closed  form  but 
must  be  determined  either  approximately  or  indirectly.  It  is  difficult  to  find  an 
approximate  form  of  q (0)  which  is  very  accurate  over  a  wide  class  of  nozzle  geometries. 
By  indirect  means  the  exact  form  of  q  versus  4>  may  be  determined  and  a  table  may 
be  constructed.  This  would  proceed  in  the  following  manner: 


4>  -  2k(0  -  0th)  =  2 K  J*qdz'  , 

where  z  is  the  axial  coordinate.  Furthermore,  for  circular  cross-sections,  the 
continuity,  isentropic,  and  isoenergetic  relations  may  be  used  to  relate  the  chamber 
radius  r  to  the  velocity  q  .  The  result  is 

^  ^7  + 1  7-1  -1y/{2(y-i)} 

This  equation  could  be  used  directly  for  given  nozzle  geometry,  to  calculate  q($)  . 
This  procedure,  however,  is  rather  clumsy  from  the  point  of  view  of  numerical  calcula¬ 
tions,  for  several  reasons  that  we  shall  not  discuss. 
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Instead,  a  much  simpler  procedure  was  employed  which  results  in  a  straightforward 
computational  procedure.  By  differentiation  of  the  above  relations,  a  first  order 
differential  equation  is  found  to  describe  Jq (4>)  .  The  equation  can  merely  be  solved 
simultaneously  with  the  equations  for  £  ,  ,  £<2>  ,  and  f3  .  The  equation  is 

of  the  form 

dq  dq  dr  dz  d 4>  1  dr 

d <t>  dr  dz  d</>  d 4>  2«q  dr/dq  dz 


so  that  it  only  remains  to  find  an  analytical  form  for  dr/dz  .  For  typical  geometries, 
this  is  a  simple  matter. 


The  major  portion  of  the  calculations  has  been  performed  for  the  nozzle  which  we 
have  already  considered,  with  the  generatrix  shaped  as  a  circular  arc  near  the  throat 
and  with  a  smooth  transition  to  a  conical  nozzle  in  the  remainder  of  the  convergent 
portion,  as  shown  in  Figure  1.  The  differential  equation  for  this  case  is  found  to  be: 


dq 

d$> 


(2R(r  -  1)  -  (r  -  l)2)* 
2*q(R  +  1  -  r)  dr/dq 

tan  01 
2*  q  dr/dq 


in  circular  region 


in  conical  region  , 


(167) 


where  all  lengths  are  nondimensionalized  with  respect  to  the  throat  radius.  Accord¬ 
ing  to  (167),  q  and  dq/d$  (but  not  d2q/d02)  are  continuous  at  the  matching  point 
(z  =  -  Rsintfj)  between  the  circular  and  conical  portions.  This  is  necessary  and 
sufficient  to  perform  the  integration. 

a 

The  initial  condition  is  q(0)  =  1  .  Note  that  </>  =  0  is  a  singular  point  since 
r  =  1  and  dr/dq  =  0  there.  One  can,  however,  find  a  regular  expansion  which  gives 
dq/d$(0)  =  {  (in  agreement  with  the  value  (dqVd0)tj,  =  1)  thereby  allowing  the  numerical 
integration  to  be  performed.  From  this  same  series  expansion,  the  term  b  appearing 
in  the  initial  conditions  for  £  may  be  found.  For  the  present  geomentry,  it  is 


2  -  y2 
6(y  +  l) 


The  relationship  between  the  nonaimensional  distance  from  the  throat  z  and  the  mean 
velocity  q  is 


-  (2R(r  -  1)  -  (r  -  l)2)*  if  0  >  z  >  -  R  sin  9l 


tan  &, 


1  -  r  +  R 


cos  6.  -  1 


if  z  <  -  R  sin  . 


The  system  of  differential  equations  (163)  through  (167)  have  been  solved  using  a 
Milne  variable-step-size  technique  by  means  of  a  Fortran  IV  rrogram  on  an  IBM  7094 
computer.  The  calculations  are  made  with  double  precision  and  a  predictor-corrector 
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scheme  is  employed  which  gives  an  accuracy  of  seven  (or.  in  other  words,  the  ratio  of 
the  error  in  the  integration  to  the  correct  value  is  of  the  order  of  10* 7).  These 
are  found  to  be  the  optimum  precision  and  accuracy  within  the  limits  of  the  operating 
range. 

The  value  y  -  1.  2  was  taken;  co  values  were  in  the  range  0  <  co  <  10  and  s^ 
values  in  the  range  0  <  s„b  <  9  .  For  each  combination  of  oj  and  s^j,  ,  the  system 
was  integrated  along  the  nozzle  length  beginning  at  the  throat  where  the  Mach  number* 
is  unity  and  extending  to  the  point  where  the  Mach  number  is  0. 05.  The  solutions 
were  determined  at  certain  specified  values  of  the  Mach  number  by  interpolation  and 
then  the  admittance  coefficients  were  calculated  at  these  specified  points  by  means 
of  (162). 

In  using  these  results,  one  would  take  the  values  of  the  admittance  coefficients 
at  that  Mach  number  which  equals  the  entrance  Mach  number  of  the  particular  nozzle  of 
Interest.  In  this  way,  as  already  indicated  at  the  beginning  of  this  section,  one 
integration  presents  information  for  an  infinity  of  contraction  ratios.  Of  course, 
we  are  considering  here  a  conical  nozzle  joined  directly  to  a  cylindrical  chamber. 

In  practice,  a  smooth  transition  would  occur  between  the  conical  and  cylindrical 
portions.  However,  an  exact  calculation  of  this  actual  situation  would  require  one 
integration  for  the  contraction  ratio.  So  the  approximation  employed  in  the  calcula¬ 
tions  is  a  large-time-saving  technique.  The  validity  of  this  approximation  was  verified 
by  calculations  with  the  more  realistic  nozzle  shape  of  Figure  17.  These  calculations 
are  discussed  in  Section  15. 

The  radius  of  curvature  of  the  circular  portion  of  the  nozzle  wall  was  equal  to 
the  throat  diameter  (R  =  2)  while  the  semi-angle  of  the  conical  portion  was  30°  in 
these  calculations  (Figure  1  is  actually  drawn  with  these  specifications).  It  is 
important  to  remcaber,  however,  that  the  admittance  coefficients  may  be  applied  by  a 
scaling  procedure  to  a  whole  class  of  nozzle  shapes.  As  shown  in  Section  13  and 
illustrated  in  Figure  2,  any  nozzle  shape  which  can  be  obtained  by  "stretching”  the 
wall  of  the  reference  nozzle  uniformly  in  the  axial  direction  by  a  factor  1/5  is  a 
member  of  this  class,  “arinking”  the  nozzle  shape  is  simply  considered  by  letting  J3 
be  larger  than  unity.  The  scaling  rules  are  given  in  the  figure.  Scaling  in  the 
ladial  direction  is  trivial,  since  all  lengths  are  nondimensionalized  with  respect  to 
the  throat  radius.  Note  that  the  radius  of  curvature  at  the  throat  R  and  the  con¬ 
vergence  angle  dx  do  not  vary  independently  in  the  scaling  procedure.  So  all 
“conical”  nozzles  cannot  be  obtained  exactly  by  scaling  from  any  one  reference  nozzle. 
However,  for  the  range  of  shapes  of  practical  Interest,  good  approximations  can  usually 
be  obtained  by  scaling  from  the  reference  nozzle.  This  is  discussed  in  detail  in 
Section  19. 

Now  that  the  possibility  of  scaling  has  been  included,  has  been  replaced  by 

svh/5  and  this  value  will  no  longer  necessarily  correspond  to  an  eigenvalue.  Therefore, 
rather  than  taking  the  usual  eigenvalues  for  svb  ,  it  is  more  convenient  to  use 
integral  values.  If  results  are  desired  for  non- integral  values  of  s„h  ,  interpolations 
may  be  made  on  certain  convenient  cross-plots. 


Note  that  the  mean  flow  Mach  number  may  replace  the  axial  coordinate  as  the  independent 
variable  for  the  purpose  of  presenting  the  results. 
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Figures  3  through  6  show  results  of  the  numerical  integration  for  a  sample  case 
(o>  =  0.5,  s„i,  =  1.0)  where  the  frequency  is  not  so  high.  In  Figure  3,  £r  is  plotted 
versus  axial  distance  showing  a  gradual  change  in  Cr  due  to  the  relatively  large 
pressure  wavelength.  Figure  4  shows  to  be  undulating*  rapidly  due  to  the 

relatively  small  entropy  and  vorticity  wavelengths.  Note,  of  course,  that  pressure 
waves  propagate  with  the  speed  of  sound,  while  entropy  and  vorticity  waves  propagate 
with  the  subsonic  gas  velocity.  Figures  5  and  6  show  the  admittance  coefficients  0^. 
and  8^  (which  are  the  most  pertinent  from  a  stability  point  of  view)  plotted  versus 
axial  distance.  Superimposed  upon  a  gradual  change  due  to  the  pressure  waves,  we  see 
a  rapid  undulation  due  to  the  entropy  and  vorticity  waves.  It  is  apparent  that  the 
admittance  coefficients,  and  therefore  the  stability  of  the  motor  operation,  can  be 
quite  sensitive  to  contraction  ratio  variations. 

At  higher  frequencies  tho  oscillations  become  more  severe,  since  undulations  in 
the  admittance  coefficients  occur  due  to  pressure  waves  in  that  case.  The  undulations 
due  to  entropy  and  vorticity  waves  become  still  more  rapid.  It  is  worthwhile  to 
discuss  the  nature  of  the  phenomenon  at  higher  frequencies  and  its  effect  upon  the 
utility  of  the  calculations. 

An  interesting  phenomenon  occurs  in  the  solution  of  the  second  order  equation  for 
$  and,  therefore,  in  £  ,  |(1)  ,  £<2)  .  and  the  admittance  coefficients.  The  real 
part  of  the  coefficient  of  $  in  the  second  order  Equation  (158)  is 

w2  -  qcz/(r'l)s2h  =  4*2j  . 

This  may  become  negative  for  low  values  of  the  ratio  u/s,vh  ,  especially  in  the  high 
Mach  nurber  region  (where  q  is  large).  Note  that  this  never  happens  for  longitudinal 
oscillations  where  =  0  •  If,  for  a  given  oscillation  (co  and  are  fixed), 
this  number  changes  sign  as  q  varies  with  0  ,  the  equation  has  a  turning  point. 

In  the  region  of  low  q  ,  where  this  number  is  positive,  the  solution  for  $(<£)  is 
undulatory,  while  in  the  region  of  high  q  where  the  number  is  negative,  $(0)  is 
not  undulatory.  Undulations  in  $(0)  imply  that  longitudinal  undulations  as  well  as 
transverse  undulations  occur  in  the  nozzle. 

This  turning  point  occurs  in  the  range  of  physical  interest  for  either  transverse 
mode  undulations  in  configurations  where  the  nozzle  volume  is  large  compared  to  the 
chamber*  volume  or  for  mixed  mode  undulations.  The  physical  reason  is  that  the  same 
frequency  may  be  a  lew  frequency  in  the  small  diameter,  high  q  ,  region  near  the 
throat  while  it  is  at  the  same  time  a  high  frequency  in  the  large  diameter,  low  q  , 
region  away  from  the  throat.  In  the  case  of  transverse  mode  undulations  in  configura¬ 
tions  where  the  nozzle  volume  is  large  compared  to  the  chamber  volume,  the  frequency 
is  determined  by  a  characteristic  dimension  which  is  smaller  than  the  chamber  radius. 
Similarly,  the  characteristic  dimension  for  mixed  mode  undulations  is  smaller  than 
the  chamber  radius.  So,  in  both  cases,  the  nozzle  may  have  two  separate  regions. 

One  region,  where  the  radius  is  large  compared  to  the  characteristic  dimension,  con¬ 
tains  longitudinal  undulations  as  well  as  transverse  undulations.  This  region 
obviously  will  always  occur  for  mixed  mode  undulations,  but  even  in  that  case  need 


*  Undulation"  pertains  to  3pace-wise  variations  while  “oscillation"  pertains  to  tiae-w<se 
variations. 

t  Aa  already  indicated,  the  chamber  is  that  portion  which  is  of  constant  dianeter  while  the 
nozzle  is  the  convergent  portion. 


not  extend  all  the  way  to  the  throat.  The  other  region,  where  the  radius  is  small 
compared  to  the  characteristic  dimension,  contains  only  transverse  undulations.  This 
region  always  occurs  for  transverse  mode  undulations  and,  in  that  case,  unless  the 
nozzle  volume  is  quite  large  compared  to  the  chamber  volume,  it  extends  over  the  whole 
volume.  It  is  seen,  therefore,  that  in  a  given  situation  either  or  both  regions  may 
occur. 

Any  undulation  in  \  means  that  $hr  ^and  #hl  pass  through  zero,  although,  in 
general,  not  simultaneously.  This  means  £  can  become  large  locally  with  the  result 
that  ,  and  the  admittance  coefficients  also  exhibit  wild  behavior  locally. 

These  regions  of  wild  behavior  do  appear  in  the  results  of  the  calculations  for  certain 
ranges  of  high  <4/s„h  and  low  Mach  number.  The  simple  physical  interpretation  is 
that  nodal  points  occur  in  these  regions  and  the  reciprocal  of  the  admittance  co¬ 
efficient  (which  is  itself  an  admittance  coefficient)  is  small.  So,  the  wild  behavior 
in  the  mathematics  is  seen  not  to  be  a  physical  wild  behavior  by  merely  observing  the 
reciprocal  of  the  results. 

However,  since  wild  behavior  means  large  changes  over  a  small  range,  it  is  not 
feasible  to  interpolate  between  two  calculated  values  in  this  range  to  obtain  a  third 
value*.  It  is  desirable  therefore  to  find  another  means  for  determining  the  admittance 
values  in  this  troublesome ^range.  This  has  been  accomplished  by  the  development  of 
asymptotic  solutions  for  l  .  ,  #<2)  .  and  f3  which  apply  to  the  low  Mach  number 

range  of  the  conical  nozzle.  This  asymptotic  theory  is  discussed  in  Section  17. 

In  addition  to  the  admittance  coefficients  G  ,  8  ,  and  C  ,  two  other  complex 
admittance  coefficients  are  useful  and  have  been  calculated.  They  are 

fy  +  g _ £ 

a  \  2  J  £*/(■/- n  +  i^/2 

e  =  ra-J-. 

la) 

c 

where  <oc  is  the  frequency  nondimensionalized,  as  in  Part  I,  by  the  ratio  of  the 
steady-state  stagnation  speed  of  sound  to  the  nozzle  entrance  radius. 

a  is  the  admittance  coefficient  to  be  used  in  the  relation  U  =  (a/y) P  at  the 
nozzle  entrance  in  the  absence  of  vorticity  and  entropy  perturbations.  This  is  essen¬ 
tially  identical  to  (105)  and  (107),  except  that  the  separation  of  variables  scheme 
and  nondimensional  scheme  are  slightly  different.  (See  the  discussion  following  (160)). 
Of  course,  when  s„h  =  0  ,  a  is  also  the  admittance  coefficient  corresponding  to  (107) 
for  isentropic  longitudinal  oscillations. 

8  is  a  combination  of  G  and  6  which  becomes  important  in  typical  combustion 
instability  applications,  as  will  be  shown  in  a  later  section  treating  the  application 
of  the  results.  It  can  be  shown  that,  for  low  Mach  numbers,  a  and  -8  are  approxi¬ 
mately  equal.  This  means  that  at  low  mean- flow  Mach  numbers  8  becomes  approximately 
independent  of  ,  even  though  G  and  8  are  dependent  upon  it. 

*  Due  to  time  and  expenses,  neither  is  it  feasible  to  perform  the  calculations  with  a  sufficiently 

fine  mesh  in  the  parameter  space  to  allow  interpolation. 
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The  results  of  the  calculations  for  the  real  and  imaginary  parts  of  the  five  co¬ 
efficients  Q,  8,  C,  £,  and  Ct  are  presented  as  functions  of  frequency  co  with 
svh  and  Mach  number  as  parameters  in  Figures  7  to  11.  Here  the  frequency  is  non- 
dimensional  ized  by  the  ratio  of  the  steady-stage  stagnation  speed  of  sound  to  the 
nozzle  throat  radius.  In  order  to  obtain  o>c  .  o>  must  be  multiplied  by  the  square 
root  of  the  area  contraction  ratio. 

One  of  the  most  interesting  results  is  that  the  nozzle  may  have  a  destabilizing 
effect  upon  the  transverse  modes  of  oscillation.  This  is  indicated  by  negative  values 
of  the  real  part  of  a  and  positive  values  of  the  real  part  of  £  .  (The  importance 
of  the  signs  of  these  quantities  will  be  discussed  in  Section  19. )  Negative  0^  and 
positive  £r  generally  seem  to  occur  in  the  range  of  "purely”  transverse  modes  where 
o)c  is  close  to  svb  .  So  here  the  nozzle  would  have  a  destabilizing  effect.  For 
longitudinal  modes  and  those  mixed  modes  where  the  longitudinal  dimensions  are  most 
significant  in  determining  the  frequency  (<ac  »  s^)  ,  0^  is  positive  and  £r  is 
negative,  so  that  the  nozzle  has  a  damping  effect  upon  the  oscillations. 

It  is  seen  from  the  figures  that  the  wild  behavior  of  the  solutions  at  low  Mach 
numbers  makes  interpolation  of  the  results  most  difficult.  For  tnis  reason,  as  already 
mentioned,  the  asymptotic  solution  has  been  employed  at  these  low  Mach  numbers  in 
order  to  determine  the  admittance  coefficients  in  that  range. 

The  figures  indicate  that  the  value  of  the  admittance  coefficient  C  are  generally 
quite  small  compared  to  the  coefficients  &  and  8  .  This  and  the  fact  that  the 
amplitude  of  the  entropy  oscillation  is  small  compared  to  the  amplitude  of  the  pressure 
and  velocity  oscillations  in  most  situations  of  physical  interest  mean  that  usually 
(159)  may  be  simplified  to  the  following  form 

U  +  GP  +  8U  =  0  .  (159a) 

As  shown  in  Figure  8,  8=0  whenever  svll  =  0  .  Furthermore  also  equals  zero 

and  it  follows  that  (-cx/y)  and  G  are  identical  in  that  case. 

It  should  be  noted  that  the  results  of  the  calculations  for  the  standard  three- 
dimensional  axisymmetric  nozzle  may  be  scaled  for  use  with  certain  annular  nozzles. 

The  major  restriction  is  that  the  inner  wall  of  the  annular  nozzle  must  have  the  same 
shape  as  a  stream  tube  contour  in  the  three-dimensional  nozzle.  This  implies  that 
the  two  nozzle  flows  are  identical  in  the  steady  state  (that  is,  of  course,  only  in 
the  common  region  where  both  flows  exist).  Also,  under  the  long-nozzle,  one-dimensional 
steady-state  flow  assumption,  this  means  that  the  ratio  of  the  outer  wall  radius  to 
the  inner  wall  radius  is  constant  along  the  convergent  section  of  the  nozzle. 

The  equations  for  the  annular  nozzle  may  be  separated  in  the  same  manner  and  the 
same  differential  equations  remain  to  be  integrated  as  in  the  three-dimensional  case. 
However,  now  s^  is  no  longer  the  hth  root  of  J'(x)  -  0  but  rather  it  is  the 
hth  root  of  J'(x)Y'(/3x)  -  J'OSx)Y'(x)  =  0  .  Here  3V  and  Xv  are  Bessel  functions 
of  the  first  and  second  kind,  respectively,  and  /3  is  the  ratio  of  the  inner  to  outer 
wall  diameter.  ( v  is  an  integer,  here. )  So  using  the  proper  value  of  s„b  ,  the 
results  of  the  three-dimensional  nozzle  calculations  for  both  admittances  and  flow 
properties  may  be  used  for  the  annular  nozzle.  The  values  of  s^  for  various  annuli 

may  be  found  in  Reference  10*. _ _ 

•  In  that  Reference,  £  Is  defined  as  the  reciprocal  of  our  £  ,  so  that  their  value  of  Sj,h 
■ust  be  Multiplied  by  their  #  to  obtain  our  value  of  s^  . 
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A  United  number  of  calculations  have  been  performed  wherein  the  throat  wall  curva¬ 
ture,  the  cone  angle,  and  the  ratio  of  specific  heats  have  been  changed.  It  was  found 
that  changing  the  last  parameter  from  the  standard  values  of  7  =  1.  2  to  y  =  1. 4 
generally  produced  a  change  in  the  admittance  coefficients  of  only  a  few  percent.  The 
other  two  parameters  affected  the  results  more  significantly,  as  shown  in  Figures  12 
through  16  where  the  results  for  different  values  of  these  parameters  are  compared. 

Calculations  were  made  with  R  =  3.0  (versus  R  =  2. 0  in  the  standard  cases)  and 
=  30°  and  also  with  R  =  2. 0  and  8^  =  15°  (versus  8X  -  30°  in  the  standard 
case?).  When  R  was  changed  and  8{  left  constant,  the  results  changed  most  signi¬ 
ficantly  in  the  high  Mach  number  range  near  the  throat.  Further  upstream  in  the  low 
Mach  number  range,  the  difference  between  the  R  =  2. 0  and  R  =  3. 0  cases  is  smaller 
On  tne  basis  of  this  small  amount  of  evidence,  it  seems  that  far  away  from  the  throat 
the  results  do  not  depend  very  strongly  on  the  particulars  of  the  nozzle  shape  near 
the  throat.  When  8 t  was  changed  and  R  left  constant,  the  solution  near  the  throat 
did  not  change,  of  course.  Only  in  the  conical  portion  of  the  nozzle  was  a  change 
produced. 


15.  EFFECT  OF  TRANSITION  REGION  BETWEEN  CYLINDRICAL 
CHAMBER  AND  CONICAL  CONVERGENT  NOZZLE 

As  mentioned  in  the  previous  section,  it  is  expedient  to  disregard  the  actual  shape 
of  the  nozzle  entrance  portion,  by  means  of  the  approximate  assumption  that  the  conical 
portion  and  the  chamber  are  directly  connected.  With  this  approximation  the  admittance 
coefficients  for  an  infinity  of  contraction  ratios  may  be  obtained  from  only  one 
numerical  integration. 

In  order  to  ascertain  whether  this  approximation  does  indeed  produce  the  negligible 
error  which  would  be  expected,  coefficients  were  calculated  for  a  very  limited  number 
of  cases  with  more  realistic  nozzle  entrance  portions.  As  Figure  17  indicates,  these 
calculations  have  two  phases.  The  first  ohase  involves  the  deteimination  of  the 
admittance  coefficients  at  the  entrance  to  the  actual  nozzle  with  contour  ABCD.  For 
the  sake  of  fair  conparison  with  the  results  of  the  calculations  for  the  approximate 
nozzle  shape  with  contour  ABE,  we  must  include  the  effect  of  the  cylindrical  portion 
ED.  So  the  second  phase  involves  the  calculation  of  the  admittance  coefficients  at 
the  E-end  of  the  cylinder  given  the  coefficients  at  the  D-end  (which  were  the  results 
of  the  first-phase  calculations).  It  is  the  results  of  these  second-phase  calculations 
(that  is.  the  coefficients  at  the  E-end  of  the  cylinder)  which  must  be  compared  to 
the  results  of  the  calculations  for  the  approximate  contour.  Obviously  the  entrance 
Mach  numbers  for  the  two  nozzles  shown  in  the  figure  are  identical. 

In  the  calculation  of  the  admittance  coefficients  for  the  actual  nozzle,  all  the 
differential  equations  remain  the  same  except  the  equation  for  the  velocity  q  which 
has  a  different  form  in  the  transition  region.  In  that  region  '.he  equation  becomes 

dq  _  1  (2R2(re  -  r)  -  (r  -  re)2)* 

dx  2*q  dr/dq  re  -  r  -  R2 

z  =  ze  +  (2R2(re  -  r)  -  (re  -  r)2)«  , 


>■ 


(169) 
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for  the  range  >  z  >  ze  .  where 

ze  =  2^  -  R2(2(i  -  cos  i9j)  -  (1  -  cos  0j)2)* 

R,(l  -  cos  6.)  +  R  ( 1 - ^-r-j  +  1  -  r 

2  1  cos  6J  e 

**  tan 

Here  the  transition  region  wall  contour  is  given  by  an  arc  of  a  circle  with  radius  of 
curvature  R2  and  the  radius  of  the  nozzle  entrance  cross-section  is  re  .  The 
velocity  is  still  governed  by  (167)  in  the  other  regions  of  the  convergent  portion  of 
the  nozzle. 

With  the  above  staple  aodifications  it  is  possible  to  find  the  admittance  coefficients 
for  this  nozzle  by  integration  of  the  equations.  Now  it  remains  to  calculate  what  the 
proper  value  of  the  admittance  coefficient  at  the  end  of  ‘.he  approximate  nozzle  would 
have  to  be.  This  value  then  must  be  compared  to  the  results  presented  in  the  previous 
section  in  order  to  determine  the  accuracy  of  using  the  approximate  nozzle  shape. 

This  second  calculation  involves  the  determination  of  the  admittance  coefficients 
at  one  end  of  a  cylinder  given  the  coefficients  at  the  other  end.  So,  an  analysis  of 
three-dimensional  oscillations  with  vorticity  waves  and  entropy  waves  in  a  cylinder 
with  uniform  mean  flow  is  involved.  This  analysis  is  essentially  a  specialization  of 
that  already  performed  for  the  variable  cross-sectional  area  nozzle  and  is  outlined 
in  the  following  discussion. 

The  linearized  equations  of  motion  (17)  through  (20)  may  be  separated  in  a  cylindrical 
coordinate  system  by  the  scheme  Indicated  by  (160)  with  the  addition  that 


w'  =  UD(z)  i  Jvh(s„hr)0'(0)ei“t 
P'  =  R(z)Jj,h(sl/hr)9(0)ei"t  . 

The  reference  length  in  the  nondimensional  scheme  is  now  taken  as  the  radius  of  the 
cylinder  walls,  z  =  0  is  one  end  of  the  cylinder  and  z  =  ■i  is  the  other  end.  A 
positive  uniform  mean  flow  q  in  the  axial  direction  exists  in  the  cylinder. 


The  resulting  system  of  ordinary  differential  equations  may  be  solved  to  obtain 
the  following: 

-  iX.z  „  iX,z 

P(z)  =  <Xe  ‘  +/3e  2 


U(z)  =  - 


a  1  iXjZ  0  iX,z 

y  0}  +  qXj  7  0}  +  qX2 


e  + 


(■=?)' 


1  + 


(>v) 


fiS35SSS»#K3l<* 


(170) 


msxm  qjym  to 
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: 


U(z)  =  >J(z)  -  i 


a  e 


iXjZ 


i\z 


■  +  -- 


y  co  +  qXj  y  a)  + 


-  iS 


ojq 


at  , 
-  i—  z ! 
»  Q 


(170) 


Cxi2  +  q^s2 


vh 


ll 

-  1-  z 

ffl(z)  -  ye  e  Q 


R(z)  =  ~  [P(z>  -  i(z)]  , 

where  a  ,  ft  ,  S  ,  and  e  are  undetermined  constants  and 

\.2  =  {*>Q  ±  b>2  ~  svh  +  -  7)  • 

In  the  case  of  interest  the  admittance  coefficients  at  z  =  0  are  known;  that  is, 
the  complex  numbers  A  ,  B  ,  and  C  are  known  in  the  relation 


U(0)  +  GP..0)  +  8U(0)  +  Ci4(0)  =  o 


(171) 


by  means  of  the  first-phase  calculation  for  the  nozzle. 

first  by  applying  (170)  to  determine  U(0)  ,  P(0)  ,  l)(0)  ,  and  £( 0 )  in  (171) 
.H  ■  *en  using  (170)  to  obtain  four  relations,  one  identity  for  each  of  U(-£)  ,  P("t)  , 
U(-t)  and  M£)  ,  we  have  a  system  of  five  linear  equations  for  the  four  unknown 
constants  a  ,  /3  ,  8  ,  and  e  This  means  a  certain  relation  must  exist  between 
the  flow  properties  at  z  =  -i  .  The  relation  is  readily  obtained  by  setting  the 
determinant  of  the  system  of  five  equations  equal  to  zero.  It  is  of  the  form 


IK-t)  +QPtf)  +SUtf)  +C Mfo  =  0  . 


Since  the  formulas  for  Q  ,  @  ,  and  C  are  rather  cumbersome,  they  are  not  presented 
here.  The  reader  may  easily  produce  them  solely_by  means  of  (170)  and  (171).  The 
formulas  allow  the  calculation  of  3  ,  S  and  C  given  y  ,  w  ,  svh  ,  ,  q  ,  (1  , 

8  ,  and  C  .  (Note  that  0L.fi.  5  ,  and  e  have  been  eliminated. )  Due  to  the  fact  that 
different  reference  lengths  have  been  used  in  the  two  phases  of  the  calculations,  the 
above  value  of  co  is  greater  than  the  u  used  in  the  calculations  of  Q  ,  8  ,  and  C 
by  a  factor  equal  to  the  chamber  radius  divided  by  the  throat  radius.  Also  since  't 
must  correspond  to  the  length  ED  in  Figure  17,  we  see  that  ■i  =  (R2/re)  tan  (6^/2)  . 


The  calculations  were  performed  for  various  values  of  o>  ,  s^  ,  Rz  ,  and  q  = 
evaluated  at  the  nozzle  entrance  (which  determines  the  contraction  ratio),  y  -  1.2, 
dy  =  30°  ,  and  R  =  2. 0  were  taken  in  the  calculations.  Table  II  presents  sample 
results  which  show  the  error  produced  in  using  the  approximate  nozzle  shape. 


The  first  complex  number  given  for  the  coefficients  is  the  value  calculated  at  the 
entranc0  of  th“  actual  nozzle  while  the  second  number  is  the  proper  effective  value  to 
be  used  at  v.he  entrance  of  the  approximate  nozzle.  (This  number  resulted  from  the 
second  phase  of  the  calculations. )  The  third  number  is  the  result  of  the  calculations 
for  the  approximate  nozzle  discussed  in  Section  14.  So  the  error  is  seen  by  comparing 
the  second  and  third  numbers.  Note  that,  instead  of  comparing  the  real  and  imaginary 
parts  separately,  one  should  more  properly  compare  magnitudes  and  phases. 
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In  the  first  four  examples  the  frequency  oj  is  large  compared  to  the  eigenvalue 
8vh  .  In  those  cases  the  frequency  is  primarily  determined  by  the  longitudinal 
dimension.  The  solutions  are  varying  slowly  through  the  transition  region  as  indicated 
by  the  small  difference  between  the  first  two  numbers.  Most  importantly,  the  small 
difference  between  the  second  and  third  numbers  shows  that  the  error  introduced  by  the 
approximate  nozzle  is  small  even  when  the  radius  of  curvature  R2  is  not-so-snall. 

In  the  next  four  cases  co  and  svh  are  equal.  Since  the  frequency  is  nondimension- 
alized  with  respect  to  the  throat  radius,  the  chamber  frequency  o>c  is  actually 
larger  than  the  eigenvalue.  This  is  the  region  of  mi~ed  mode  oscillations  where,  for 
the  flow  velocities  indicated,  the  radial  and  longitudinal  dimensions  are  equally 
important  in  determining  the  frequency.  Here  the  flow  properties  may  vary  more  rapidly 
through  the  transition  region,  resulting  in  the  possibility  of  significant  errors  by 
use  of  the  approximate  nozzle.  Still,  as  shown  in  Table  II,  the  error  generally  is 
satisfactorily  small  for  small  radii  of  curvature  at  the  nozzle  entrance.  This  is 
expected,  since  the  actual  nozzle  approaches  the  approximate  nozzle  as  R2  tends  to 
zero. 

The  last  example,  when  compared  with  the  one  immediately  preceding  it.  indicates 
the  improvement  in  the  accuracy  of  the  approximation  as  the  frequency  decreases  below 
the  eigenvalue.  There,  the  range  of  “pure"  transverse  oscillations  is  approached. 

In  conclusion,  the  calculations  show  that  the  approximation  usually  introduces 
little  error.  However,  in  those  ranges  where  the  flow  properties  vary  rapidly  in  the 
axial  direction,  significant  errors  may  be  obtained  for  large  radii  of  curvature  at 
the  nozzle  entrance. 


16.  FLOW  PROPERTIES 

It  is  interesting  to  examine  the  actual  velocities,  pressure,  and  entropy  of  the 
oscillatory  nozzle  flow,  even  though  the  knowledge  of  these  quantities  is  not  required 
for  the  determination  of  the  admittance  coefficients.  The  determination  of  the  flow 
properties  should  lead  to  a  better  physical  understanding  of  the  nature  of  the  oscilla¬ 
tion  in  the  nozzle. 

In  order  to  determine  these  quantities  one  must  solve  (158)  for  $  and  substitute 
the  solution  into  (102).  In  addition  f0  ,  ft  ,  and  f2  must  be  determined  and  cr  . 

Cj  .  and  C2  must  be  specified.  In  order  to^solve  (158)  it  was  convenient  to  first 
change  the  independent  variable  from  4>  to  4>  =  2 x(<p  -  <j>tb)  .  The  differential 
equation  becomes 


d2$  d$ 

b  7T  "  +  ih)  ~y+  (j  -  ik)$  =  -  (G  +  H)  ,  (172) 

d <p  d  <P 


where  the  coefficients  b  ,  g  .  h  ,  j  ,  and  k  are  those  defined  immediately  following 
(165)  and  G  and  H  are  given  by 


G 

H 


(0-/2*)  erf. 


df,  s*bc2/(r'l)" 


2  io) 


1  ♦ 
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It  is  convenient  for  the  purpose  of  nuaerical  integration  to  reduce  this  second-order 
complex  differential  equation  to  four  first  order  equations  by  means  of  the  following 
definitions:  $  =  y2  +  iy2  .  d§/d0  =  y3  +  iy„  .  G  =  Gr  +  iGj  .  and  H  =  Hr  +  iHj  . 

Then  the  system  of  equations  becomes 


y'  =  [gy3  -  jy l  -  hy,  -  ky2  -  Gr  -  Hr]/b 


(173) 


yj  =  [gy„  -  jy2  +  hy3  +  kyt  -  Gj  -  Hj]/b  . 


/ 


Since  b  =  0  at  the  initial  point  of  the  integration  (as  before),  one  must  insure 
that  the  regular  solution  is  obtained  by  expanding  the  solution  in  a  Taylor  series 
about  the  singular  point  at  the  throat  in  order  to  determine  the  value  of  the  solution 
and  its  first  derivative  at  that  point.  These  values  are  then  used  in  the  first  step 
of  the  numerical  integration.  Again,  we  use  the  steady-state  velocity  profile 
q2  =  1  +  0  +  b 4?  +  .  .  .  So.  by  expanding  the  coefficients  and  inhoaogeneous  parts 
of  (173)  in  a  series  in  <£  and  determining  the  coefficients  in  the  series  solution 
for  (173)  in  the  vicinity  of  the  throat,  one  finds  the  following  initial  conditions*: 


y3(0) 

=  i  , 

y{(0) 

=  Br 

y2(0) 

=  o  , 

y'(0) 

=  Bi 

y3<°> 

=  Br  • 

*3<°> 

=  *, 

y,(Q) 

=  B,  . 

y'(0) 

=  2C, 

where  the  following  definitions  apply 


(174) 


In  linearized  problems  dealing  with  periodic  phenomenon  the  amplitude  and  time-phase  are 
arbitrary  in  tbe  final  solution.  So  with  no  loss  of  generality  4>(0)  =  1  ras  chosen. 
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2<£iu  -  2(r  +  1)1, 
4p  +  (y  +  1)2] 


C,  = 


2wl,  +  2(7  +  1)I„ 


1  4  [o>2  +  (yf  1)2] 


Defining  further,  we  have 


/VO  aO 

_  Kh  ~  ®*n 


cxr  " 


2x(y  +  l) 


y  -  i  *  s, 


co  +  -&C,t  - 


4  11  2*  (y  +  1) 


7  +  1 


+  (7  +  l)b 


Br  + 


^2 

®vh  -  w 


S*h 


B_  -  (y  -  1) 

4  r  8 


0CU  -  (7  -  l)Clr]  + 

y2  -  67  -  3  cr.  ary  cr. 

+ - - - c.  + - *■ 

2(7  +1)  2k  7  +  1  4k 


lu  =  ~ 


7  +  1 


+  (7  +  l)b 


Bi  + 


CO  —  8 


vh 


Bi  - 


3  +  7  ~  72  -  1  -  7  -  1  «" 

- wB - a) - «b  - 


y2  -  67  -  3  cri  coy  ar 


~  (7  -  l)Ctl  -  wClr  -  — - - n - =-  +• 

8  lr  (7  +  1)  4  7  +  1  4k 


When  solving  the  system  (173)  and  (174),  f0  ,  tl  ,  f2  ,  and  f3  are  determined 
simultaneously.  They  are  governed  by  the  following  system  of  equations  which  may  be 
obtained  from  (77),  (81),  (83\  and  (85): 


*ir«» 


<•>  f0l(4>) 
2  q* 


«  for«» 


fi‘(«  = 


for(°>  =  1 


foi(0)  =  0 


(175) 
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1  dq5 

-  f  nr 


flr(0)  =  0  , 


f:,(0)  = 


1  dq5 
-  f„,  — r- 


fn<(0)  =  0 


f'r(0) 


f^(0)  = 


(C^f3);  =  - 


«f  21(0) 

flr 

2  q2 

2*q2 

^2r<*>  + 

fii  . 

•2q2 

2*? 

^f3i 

1  dq2 

2^ 

2  d<£ 

f2r(0)  =  0 


f,,(0)  =  0 


czf,_(0)  =  0 


<?*•>{ 


^  ^3r  . 
2q?  * 


?f3l(0)  =  0 


The  numerical  integration  of  Equations  (175)  is  simpler  than  the  numerical  evaluation 
of  tho  integrals  appearing  in  their  exact  analytical  solution.  The  flow  properties 
were  determined  for  the  30°  conical  nozzle  described  in  Section  14.  So  (167)  and  (168) 
were  solved  simultaneously  with  Equations  (175).  The  integrations  were  performed  with 
an  IBM  7094  computer  in  similar  fashion  to  those  previously  mentioned. 

Substituting  for  $  and  d$/d0  =  2*d$/d0  in  (102),  we  obtain,  after  certain 
manipulations, 

Umod  =  2***3  +  ] 


Uarg  =  arc  tan  (yu/y3) 


Pmod  =  [*^2  +crrfir  "  aifii  “  2*Q2y3)2  + 

li 

+  (-  2*q\,  -  *<2yx  +  aTf  u  +  o-jf  lr)2Jl 

^rfn  +  aif-.r  -  *"y,  -  2*7yK 


p  =  tan 5£ll -  ?..jb 

arg  crrflr  -  Ojfjj  +  KU y2  -  2*q2y3 


V  —  w 
mod  mod 


r— 

[*yi  “  ^lr^cr  +  Cti*oi  “  °r^2r  +  °i^2p  + 

-.1 

+  (y2  -  Clrf0£  -  Cjjfor  -  crrf2i  -  cri^2r'^  J 


V  =  W 
arg  arg 


arc  tan 


*2  ~  cir^oi  ~  ^ll^or  _  ^r^l  _  cr1^2r 
*1  "  Cir*or  +  Ciifoi  ”  CTrf2r  +crif2i 


(  (176) 


,  -WTMMWggMw vwm— 
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mod 


(°r  +  * 


arg 


=  arc  tan 


Vol  +  alfor 
Vor  "  CTifoi 


where,  for  example,  P  = 


P«od  <lparg> 


(176) 


The  results  given  in  (176)  are  nondimensionalized  with  respect  to  static  throat 
conditions  and  are  consistent  with  the  variable  separation  scheme  given  by  (44).  It 
was  desired  to  have  the  final  results  nondimensionalized  with  respect  to  stagnation 
conditions  and  to  be  consistent  with  the  variable  separation  scheme  given  by  (160). 
This  was  readily  achieved  by  multiplying,  in  (176), 


Pmod 


2 y  P  y  -  l  — ' "li/(7-i) 


7 


tL.  i -i— i? 

+  1  [_  y+l 


2  V 

by  - - -  q 


y  +  v 


and 


vB0d  “W*  wmnH  by 


mod 


2  \* 


7  +  I J 


As  already  specified,  r  in  (160)  (as  well  as  in  the  equation  for  w'  and  p'  given 
in  Section  15)  must  be  interpreted  as  the  radius  divided  by  the  local  wall  radius. 


In  order  to  readily  compare  the  results  for  different  cases,  a  phase  was  added  to 
the  arguments  of  (176)  such  that  Parg(0)  always  was  equal  to  zero.  Therefore  the 
phase  angle  A  to  be  added  was  given  by 


A 


-  arc  tan 


-(2Bj  +  w) 


-  2B„ 


Sample  results  of  the  flow  properties  calculations  are  given  in  Figures  18,  19, 
and  20  for  the  case  of  an  oscillations  with  s„h  =  1.84129  and  co  =  l.  00  .  Here  the 
transverse  variation  corresponds  to  that  of  the  first  tangential  mode.  The  same 
nozzle  shape  as  discussed  in  Section  14  was  considered.  Both  cr  and  Ct  are  real 
and  equal  to  0.5  for  this  case. 


In  Figure  18,  it  is  seen  that  there  is  a  gradual  axial  variation  of  the  pressure 
due  to  the  acoustic  wave  phenomenon  in  the  nozzle  but  there  is  essentially  no  variation 
due  to  the  presence  of  entropy  and  vorticity  waves.  The  axial  component  of  velocity 
has  a  gradual  variation  due  to  the  acoustic  wave  plus  a  more  rapid  undulation  of  smaller 
amplitude  due  to  the  entropy  and  vorticity  waves,  as  shown  in  Figure  19.  The  trans¬ 
verse  components  of  velocity  have  rather  severe  undulations,  due  to  the  entropy  and 
vorticity  waves,  as  indicated  by  the  ’'ariation  of  the  radial  velocity  component  shown 
in  Figure  20. 


These  entropy  and  vorticity  waves  travel  at  the  particle  velocity,  which  is  subsonic, 
so  that  their  wavelength  is  shorter  than  the  acoustic  wavelength  but  increases  as  the 
wave  moves  towards  the  throat.  The  values  for  cr  and  Cj  in  the  figures  were  taken 
to  be  considerably  higher  than  those  normally  found  in  practice  for  the  purpose  of 
demonstration.  The  results  are  qualitatively  similar  for  a  wide  range  of  values  of 
these  complex  numbers.  Even  though  the  pressure  does  not  exhibit  the  presence  of 
entropy  waves,  the  density  and  temperature  would. 

Note  that  the  admittance  coefficients  are  independent  of  the  values  of  i /  and  C1  , 
and  therefore  independent  of  the  intensity  of  the  entropy  and  vorticity  waves.  Of 
course,  as  cr  and  C,  go  to  zero,  the  general  admittance  relation  can  be  shown  to 
reduce  to  the  irrotatlonal  admittance  relation. 


17.  ASYMPTOTIC  BEHAVIOR  OF  THE  ADMITTANCE  COEFFICIENTS 

Due  to  the  wild  behavior  observed  in  the  numerical  solutions  of  Equations  (163), 
(164),  (165),  and  (166)  in  the  low  Mach  number  range  of  the  nozzle,  it  was  desirable 
to  obtain  some  analytical  prediction  of  the  solutions  to  those  differential  equations. 
This  was  achieved  by  developing  solutions  which  were  asymptotic  in  the  sense  that 
they  apply  in  the  limit  as  Mach  number  goes  to  zero.  These  solutions  are  satisfactory 
for  the  present  purpose  since  it  is  in  the  low  Mach  number  region  (away  from  the 
throat  location)  that  the  behavior  must  be  examined. 

In  the  low  Mach  number  region  of  the  nozzle,  we  can  say  (using  the  continuity 
relation  and  neglecting  terms  of  order  q5)  that 

(2/(y  +  i))i/{*<r-»>)  =  r(q)»  , 


where  the  same  nondimensional  scheme  as  discussed  in  Section  14  is  employed  here. 
Noting  that  for  the  conical  portion  of  the  nozzle  dr  =  -  dz  tan  9  =  -  dqktan  £)/2*q  , 
we  find,  after  differentiation  of  the  above  relation,  that 

(y  +  A  ^ 2(y- 1>)  tan  0.  o  1 
\  2  J  2k  2q» 


This  may  readily  be  integrated  to  obtain 

q  =  (aip  +  c)2  . 

where  c  is  determined  by  specifying*  q  at  some  value  of  4>  and 

+  jy/Wi)}  tan  g, 


(177) 


a  = 


2k 


It  is  convenient  to  begin  the  development  of  the  asymptotic  solutions  by  a  trans¬ 
formation  of  the  independent  variable  which  is  suggested  by  the  form  (177). 


•  For  the  purposes  here,  specification  of  the  value  of  c  is  unnecessary. 
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Setting  y  =  a0  +  c  ,  it  readily  follows  that 


q  =  y2 


c  = 


7+1  7-1 

~~2  ~ 

d 


d  <p  a  dy 


(178) 


C  will  be  determined  by  finding  the  homogeneous  solution  to  (158)  and  then  taking 
the  derivative  of  the  logarithm  of  that  solution.  Using  (178)  to  substitute  into 
(172),  we  have  the  equation  with  complex  coefficients 


$£'  +  (A  -  iC)$Z  +  (B  -  iD)$h  =  0  . 


(179) 


where 

A 


-  8y3/ (7  +  1)  +  0(y7) 


B  = 


a>- 


1 

— -  + 


:i „  (y  *  iY/<’"‘>  (' 


2(7  +  l)a2  y“  2(7  +  l)a2  4a2  V  2 


\y*  7+1 


7  y2)  +  o(y6) 


c  = 


2co  2wy  .  8. 

+  — - —  +  0(y8) 


(7  +  l)a  (7  +  l)a 


4(7  -  1) 

D  =  ~t, — it?  y  +  o(y  ) 

a2  (7  +  1) 


The  linear  equation  (179)  with  variable  coefficients  cannot  be  solved  exactly; 
however,  there  is  the  well-known  WKBJ  method  for  solving  a  second  order  linear 
equation.  The  important  criterion  is  that  the  solutions  vary  much  more  rapidly  than 
the  coefficients  in  the  differential  equation.  This  is  satisfied  if  <£/ a  is  large, 
which  implies  that  the  axial  gradients  in  the  unsteady  state  are  larger  than  those  in 
the  -_t':ady  str.ts,  since  a  is  proportional  to  tan  9  ^  . 

The  *XBJ  method  involves  the  determination  of  the  two  solutions  in  the  form 

$h  =  ^(y)  .  (180) 


For  equations  with  constant  coefficients  the  derivative  </>'  is  a  constant.  In  the 
present  situation,  where  the  solution  varies  much  more  rapidly  than  the  coefficients, 
this  quantity  is  “nearly”  constant,  so  that  </»"  is  of  higher  order.  Therefore,  in 
the  first  approximation  we  may  neglect  so  that,  substituting  (180)  into  (179), 
we  have  ,  .  j, 

,f  A  iC  (  C2V  5 

xp<  =  —  +  —  ±  i  (b  +  — )  +  0(y5)  . 

2  2  V  4/ 
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This  first  approximation  is  differentiated  to  obtain 


r  ♦ — 2"  -itl _ " _ r  lii.i  fa* 

a(2(y+l)}i  y3  a{2(7+l)}«  7+1  4\a>, 


r+3  i/£>  /r+  i\2r/(r*l) 


+  0(y2) 


Now,  in  the  second  approximation,  the  </'"  is  kept  when  (180)  is  substituted  into 
(179),  with  the  following  result: 

A  iC  i 

</»'  = - +  —  ±  i (B  +  'p"  +  C2/4)*  ,  (182) 

2  2 

where  is  given  by  (181). 

For  the  third  approximation,  the  result  of  the  second  approximation  is  differentiated 
and  substituted  into  (182),  which  then  is  solved  for  </»'  .  Of  course,  this  iteration 
may  be  continued  indefinitely.  The  result  of  the  fifth  approximation  is  conveniently 
expressed  as  follows  (setting  <//  =  </>r  +  i</».): 

’/'r  =  -  +  K:y  +  K3y3  +  0(y5) 

y 


&  1  l  s2  (y  +  i\^1)/{^(r-i)} 


, u  -  ±  ~  *  *  ‘’va  I '  1 

1  |a{2(r  +  1)}«  y2  2  wa  \  2 


±  -  +  Kzy2  +  0(y4) 

a  (7  +  1)  2 


wh^re  the  following  definitions  apply 

■,  ■  ;(%)'? tT“ 


1  7  +  3  w  3  as2h  .  (y  + 

*2  =  2  7+1  a{2(7  +  1)}*  +  \  {2(7  *  1)}'  (“ 


i  (y  + 


+  9-  (7  +  l){2(y  +  1)>*  (i)  -  -  — yTl  (— 

8  \w/  8  aco3{2(7  +  l)}»  \  2 


3  as2h{2(y  +  l)}*  (y  +  l\r /(r_1)  3  {2(7  +  l)}‘(7  +  l'/a3 


7  +  l\2’,/('>''1) 


7  +  3  „a2s 2h  /7+  i\(2r’l)/(r-1)  ,  x7  a2 

“  T-T-T  “  11  “s?  HH  -10(7+1)2^j 

7+1  or  \  2  j  or 

i  Kh  (y+  iVy/(y-1? 

8  a>4  \  2  / 
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(183)  shows  that  two  solutions  exist  for  ,  which  means  that  <//  ,  and  therefore 
$b  ,  also  have  two  solutions. 

Integrating  (183),  we  find  that 

*pT  =  b  +  log  y  +  +  o(y6) 

2  4  •  (184) 

</»i  =  i -  c  or  -  +  d  , 

4 

where  b  ,  c  ,  and  d  are  constants  of  integration  and 
.  _ a  i 

1  ”  ft{2(y  +  1)}^  y  ~  2  ui  \  2  )  ?  + 

a(y  +l)  3 


-)  =  _ 


1  1  s^h /V+  lVy+1)/{2(y*1)} 


a{2(y  +  1)}t  y  2aa  \  2 


«y  K2y3  , 

— - -  +  -*—  +  0(y5) 

a(y  +1)  3 


It  can  be  shown  that  </^+^  and  '/'j correspond  to  waves  travelling  in  opposite 
directions;  corresponds  to  a  wave  moving  towards  the  nozzle  entrance  and 

corresponds  to  a  wave  moving  towards  the  nozzle  throat. 

Substitution  of  (184)  into  (180)  yields  the  following  results: 

=  Ae'r  exp  {i(^+)  -  c)>  +Be*r  exp  {- iCAj_)  -  d)}  ,  (18S) 

where  J  ,  5  ,  c  and  d  remain  unknown.  Since  exp  (b)  may  be  included  in  A  and 
6  with  no  loss  of  generality  we  will  consider  b  =  0  in  (184). 


Noting  that 


(186)  may  be  used  to  deduce  that 


c  ■  m 


$h/  \  d0  /  A  dy  / 


l  _  a././  5)  sin  -  (c  +  d)) 

r  “  (A/B)2  +  1  V  2(A/B)  COS  +  <//<">  -  (c  +  d)) 


{(*/§)»  i) 

-  =  .(»;♦>■  -0j-i>) _ 2  W6) 

1  2  +  (A/B)2  +  1  +  2(A/B)  cos  ('p[*)  +  'p[')  -  (c  +  d') 
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It  is  seen^that  the  unknown  constants  appear  Jn  two  combinations,  (A/B)  and  (c  +  d) 
If  £r  sad  were  given  at  some  initial  value  of  y  ,  the  constants  (A/S)  and 
(c  +  d)  could  be  calculated,  since  0r  ,  and  i and  their  derivatives  are 

known  as  functions  of  y  .  Once  these  constants  are  known,  (187)  coula  be  used  to 
determine  £r  and  at  other  values  of  y  . 

If  the  relations  given  by  (177)  and  (178)  are  used  to  transform  (166),  the  result 
is  that 


f  <?f ,)  =  7^  (?t,)  ♦  ■ 

dy  3  2  ay 4 

Upon  integration  we  have  the  exact  solution  as  follows: 

?f3  =  k  exp  (-iw/6ay3)  +  2  exp  (-io>/6ay3)  JJy^  exp  (-i^/Gay3)  dy  , 

where  y  is  a  dummy  variable. 

The  definition  is  made  that 

In(y)  =  exp  (~i<£/6ay3)  Jjp*  exp  (-ia/6ap)  dy  .  (188) 

Then  ve  have  the  following  form  for  the  solution: 

?f,  =  k  exp  (~io/6ay3)  +  2I3  ,  (189) 

where  k  is  a  constant  which  is  determined  when  f 3  is  known  at  a  certain  value  of  y  . 

It  remains  to  determine  In(y)  ;  an  asymptotic  series  may  be  developed  by  successive 
integration  by  parts.  First,  the  following  transformations  are  made 

y  = 

This  means  that  (188)  becomes 

ln  = 

Integration  by  parts  generates  the  asymptotic  series  approximation 


V  = 


n  +  4 


d)\(n+1)/3  e~lt 


dr 


( 188a) 


(190) 
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=  V  .  S+i  =  ~  +  2m) (tj  +  2a  -  1)  , 

bj  =  -  7)(r/  +  1)  .  bB+l  =  -  b>(7j  +  2m  +  1) (t?  +  2m)  . 

The  integral  which  would  remain  after  2M  +  1  integrations  represents  the  error  in 
using  the  series  approximation  in  place  of  the  exact  function  In  .  This  error  can  be 
shown  to  be  bounded  above  by  the  quantity 

(V  +  2M  -  1)(t?  +  2M  -  2) - (7?  +  1) 

E  =  (3)V>+2“ 

For  a  given  value  of  the  argument  t  ,  therefore,  an  optimum  value  of  the  integer  M 
exists  which  minimizes  the  upper  bound  E  .  This  optimum  situation  occurs  when  M 
equals  one  plus  the  greatest  integer  in  the  quantity  (t  -  t?)/2  .  The  smaller  the 
Mach  number  is,  the  larger  t  is,  and  the  larger  the  optimum  value  for  M  is.  In 
the  limit  of  zero  Mach  number,  t  is  infinite  and  the  optimum  value  of  M  becomes 
infinite. 

In  is  not  a  rapidly  undulating  function  of  y  for  the  values  of  a>  in  the  range 
of  interest.  In  particular,  I3  is  proportional  to  y7  in  its  leading  term,  so  that 
c*f3  is  essentially  given  by  its  homogeneous  solution.  That  solution  portrays  the 
effect  of  properties  which  propagate  with  the  gas  velocity  rather  than  the  sound 
velocity  since  exp  (-i«o/6ay3)  is  an  approximation  to  exp  (-io)/(l/2q5)dx)  . 

Asymptotic  solutions  may  also  be  found  for  Equations  (164)  and  (165).  Equations 
(117)  and  (178)  may  be  used  to  transform  those  differential  equations  to  the  following 
forms*: 


—  M(i>  +  £  -  (—  +  - - - ^ 

dy  ’  (7  +  1)(1  -  y 7 


M<1)  =  (l-y^f^ 
M<2>  =  (l-yu)f<2> 


(i,  ...  fz  +  iYa-y)/(y-i)  s3h  /  7-i 


7+1 


l/(7-l) 


- a  — -  + 

7+1  dy  2ico 


/ 

'y  +  l 

2iw' 

V  2 

.  2y2 


Equations  (189)  and  (190)  may  be  employed  to  simplify  these  equations  as  follows: 


•  F(l)  is  not  exactly  the  same  quantity  here  as  it  is  in  (88). 
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For  our  purposes,  it  will  be  proper  and  consistent  to  take  exp  (0(iy))  as  suffici- 
ently-slowly  varying  to  be  replaced  by  unity-  We  see  then  that  integrals  of  a  certain 
type  sun  to  produce  the  integral  in  (192).  This  type  of  integr.d  is 


r  y 


yn  exp 


•i  0 


(i<£/a) 


\6r  {2(7  +  l)}»y/ 


dy 


In  particular,  for  our  desired  accuracy  we  are  interested  only  in  the  cases  where 
n  =  1  and  n  =  -1  .  Note  that  the  exp  (0(i/y))  tern  represents  wave  propagation 
effects  and  varies  much  aore  slowly  than  the  exp  (0(i/y3))  tern,  whicn  represents 
particle  propagation  effects.  So  there  is  a  preferred  way  in  which  to  integrate  by 
parts.  One  finds  that,  to  the  desired  accuracy. 


y'1  exp 


(i<*/a)  *  r. 


\6y3  {2(7  +  l)}*y ) 


dy 


=  exp 


(i<2/a) 


V6y3  {2(7  +  i)}«y/ 


2  \* 


!-w±7n 


T  exp 


1 


(icV  a)  —  ± 


,6y3  “  {2(7  +  l)}iy ) 


dy 


=  exp 


j(i<^a)  ±  ^  +  1i 


,(y)  . 


where  (190)  gives  Ij(y)  and  I.j(y)  .  Now  (192)  yields 


M(1>  =  (l-y')f(1>  = 


\  -i<S/6ay' 

_ 


fy  +  ^  y3  +  ((X/b)  e-<c+d>  e'ia/y  +  eia/y) 

\  2  8 


M' 


.1— 1  ^ 

4y 


I 


.  ,  iV  V2  +  + 

1  +  jRj  y 


/  2  V  f(VB)  e-Uc^)g-?la/y 

+  \7  +  l)  11  W)e"1(c+d)e-*ia/y  +  1 J 


(193) 


where  the  quantity  a  has  been  defined  as 


CO 

a{2  (7  +  l)}i 
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In  developing  the  asyaptotic  form  of  (192)  for  £t2)  ,  integrals  appear  of  the  saae 
form  as  appeared  in  the  relation  for  .  In  addition  integrals  of  the  form 


5h«x)  = 


y°  exp 


jo 


iOJ 


a(2(y  +  1)}»  y 


dy 


appear.  It  can  be  shorn  that  exp  [+  iaVa{2(y  +  l)}*y]  can  be  put  into  a  form 
similar  to  (188a)  with  t  =  +  <S/a{2 (y  +  l)}*y  aj.d  i?  =  n  +  2  .  Then  it  follows  that 


*n<y>  =  3 


L/  »y| 

n-f  1 

ico 

V  +  l) 

j  exp 

a(2(y  +  l)Ky 

<y) 


where  I3n+2(J)  is  evaluated  from  (190)  by  using  the  above-mentioned  relations  for 
t  and  t)  .  The  “wiggly"  sign  over  I  has  been  used  to  denote  the  change  in  the 
argument  as  a  function  of  n  and  y  from  the  previous  usage  even  though  it  remains 
as  the  same  function  of  t  and  t)  .  Note  that  §_3  is  special  in  that  it  is  improper: 
its  amplitude  goes  to  infinity  as  1/y  .  as  y  goes  to  zero.  This  really  occurs 
only  since  we  chose  to  set  the  lower  limit  of  the  intetral  at  y  =  0  .  If  we  set  the 
lower  limit  at  some  finite  number,  say  e  .  the  indefinite  integral  evaluated  at  the 
lower  limit  is  merely  a  constant  which  may  be  incorporated  with  the  constant  which 
appears  before  the  homogeneous  solution  for  X2  .  This  Is  done  with  no  ;.oss  of 
generality  so  that  we  nay  take 

„  2a2  (y  +  1) 

V-3  *9 


\  a{2(y 


i «  n 

7TWzyJ 


exp 


ico 


a(2(y  +  1)}« 


2 
«  y 


Note  that  the  above  asyaptotic  solution  which  contains  only  two  terms  in  the  series 
is  in  fact  an  exact  solution  for  the  integral. 

It  now  follows  that 

M(2)  =  (1  -  y*)#(2> 

2<wk 


(7+  l)2 


exp  (-iav'eay3) 


2i£k 


(A/B)  e'1  (c*d>  e' 2iB/y  -  1  (7  +  if  m2 ***  (_iW6ay  5 


(a/E)  e-<c*d>e- ! 


y  +  is,  y3 


\2  exp  C-icVeay3) 

+  (y  +tKj  y3)((A/B)  e-» +  elB/y)  + 

y  +  ?k5  y3 


exp  (-io>/6ay3) 


(194) 
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I 


(X/5)  e'1^c+d>  e"2im/y  i<;>  +  I_<t> 

(A/B)  ei<c+dV*lB^  Vl 

<5k  2K3  +  K2  +  g^jcKj 

(7  +  l)2  4  \  2  / 


54  iexp  (-i«/6ay3) 

+ - ; - - -  x 

(r  +  i>y 


(194) 


(A/B)e-i(c+d>e-2lm/y  I<”>  +  I<+i 
(A/B)e'(c+d)  e-2l2/y  +  1 


y  + 


X(2-r)/(r-0  ^ 

)  4wy  X 


) 


The  plus  or  minus  in  the  superscr.vt  of  the  I  function  indicates  which  sign  is  taken 
for  t  in  (190). 


1 

\ 

i 

t 

t 

5 
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The  solutions  given  by  the  asymptotic  theory  provide  a  great  amount  of  insight  into 
the  oscillatory  flow  in  nozzles.  (186)  shows  that  #  is  the  sum  of  two  solutions, 
one  corresponding  to  upstream  wave  propagation  and  the  other  to  downstream  wa;  . 
propagation.  The  equations  following  (135)  show  that  goes  as  1/y  tv.  leading 
order.  Noting  that  d<£  =  2*rio b  ,  dy  =  ad<£  ,  and  d<£/dz  -  q  =  y2  we  find  1  'y  is 
linear  in  z  ,  the  axial  coordinate,  so  that  the  solutions  are  undulatory  with  constant 
wavelength  but  their  amplitudes  vary*.  The  radical  behavior  in  t  is  now  easy  to 
explain.  It  is  seen  that  $  behaves  ao  a  nearly- sinusoidal  function,  so  that 
l  -  ( 1/#'/ (d$/d<£)  behaves  nearly  as  a  cotangent  which  becomes  quite  large  in  certain 
ranges.  The  appearance  of  the  cosine  function  in  the  denominators  of  (187)  shows  that 
periodically  t  may  become  large. 

a  A 

It  is  especially  interesting  that  to  leading  order  £  does  not  depend  upon  sl/fa 
but  only  upon  a>  .  The  wavelength  can  be  shown  to  be  approximately  that  of  a  planar 
wave  generated  by  a  source  of  frequency  a>  .  In  a  cylindrical  chamber,  one  would  have 
found  a  larger  wavelength  since  (&>2  -  s2j,)‘  is  the  effective  factor  rather  than  co  . 

Equations  (193)  and  (194)  shew  that  and  ^(2>  are  affected  not  only  by 

acoustic  oscillations  but  by  oscillations  of  entropy  and  vorticity  as  well.  The 
acoustic  behavior  appears  through  the  exp  <±  im/y)  terms,  while  the  entropy  and 
vorticity  behavior  appears  through  the  exp  (-ia>/6ay3)  terms.  Both  £(1)  .d  £  2) 
have  sinusoidal  and  cosinusoidal  terms  appearing  in  denominators  of  their  expressions 
(193)  and  (194)  so  they  become  large  periodically  for  the  same  reason  that  £  becomes 
large. 


18.  ASYMPTOTIC  DEVELOPMENT  OF  THE  FLOW  PROPERTIES 

In  addition  to  the  development  of  the  asymptotic  theory  for  the  admittance  co¬ 
efficients,  it  is  convenient  to  develop  an  asymptotic  theory  to  determine  the  pressure, 


0  The  variation  is,  of  course,  due  to  convergence  of  the  nozzle. 
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density,  and  velocity  undulations.  This  theory  describes  the  effects  of  the  acoustic 
waves  and  the  entropy  and  vorticity  waves  upon  the  oscillatory  flow  properties,  at 
least  for  the  low  Mach  number  region  of  the  nozzle.  Of  course,  it  provides  information 
that  could  be  found  qualitatively  in  the  results  of  Section  16  only  by  means  of  hind¬ 
sight.  Quantitatively  it  gives  the  orders  of  magnitude  of  various  interesting  effects. 


The  asymptotic  solution  of  (172)  for  low  Mach  number  is  required.  The  relations 
(177)  and  (178)  are  employed  again  and  the  coefficients*  of  (172)  may  be  represented 
by  their  asymptotes  as  follows: 


b 

g 

h 

j 

k 


2  ( - - y7  +  0(yn)  . 

v  +  V 


— — —  y6 1  +  0(y10) 
7+1  ) 


(195) 


The  transformations  discussed  in  Section  14  are  used  in  conjunction  with  (77), 
(177)  and  (178)  to  yield 

f0  =  \3  exp  (i<£/6ay3) 

where 


\ 


3 


exp  (-i<3/6ay3')  exp 


(196) 


y  and  correspond  to  the  same  axial  position  somewhere  in  the  low  Mach  number 
region  of  the  nozzle.  Now.  using  the  definitions  following  (172),  together  with  (177), 
(178),  and  (196),  we  find  that  the  inhomogeneous  part  of  (172)  has  the  following 
asymptote: 


•  Actually,  these  coefficients  are  defined  immediately  after  (165). 
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Inspecting  the  coefficients  <195)  and  the  inhomogeneous  part  (197),  we  expect  to 
find  a  particular  solution*  $>p  of  the  form 


§  =  gi^/Ssy3  £  vn  £  . 

v  n=-o  “  n=  -as  “ 


(198) 


Substitution  of  (195),  (197),  and  (198)  into  (172)  allows  the  determination  of  the 
coefficients  ^  and  Bn  .  The  particular  solution  is  found  to  be 
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+  0(y7)  +  0(y10  e1^6*37  ). 

(199) 


The  solution  is  given  to  higher  order  in  the  undulating  terms  than  the  others,  since 
their  gradients  are  large  and  consistent  accuracy  in  the  gradients  is  required.  The 
above  solution  when  added  to  C2$h  ,  where  is  given  by  (186),  represents  the 
asymptotic  solution  of  (172),  so  we  write 


$ 


CA  +  %  ■ 


Use  of  (81),  (83),  (177),  (178),  (196),  and  the  transformations  discussed  in 
Section  14  yields 

fj  =  K  *  2\3  eii/6&y3  In (y)  =  +  0(y7  e ii/6ay3) 

f  =  X  ei"/6ay3  _22l\ 


where 


COK 


1  I  d?  1 

2  f°  "d0"d 


(200) 
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•  Clearly,  th(|  particular  solution  must  essentially  equal  the  particular  solution 

+  &^2'  of  (88);  the  only  difference  is  in  the  nondimensionalization  scheme,  as 
explained  at  the  beginning  of  Section  14. 
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The  relations  (102)  give  the  formulas  for  the  axial  dependence  of  the  flow  properties 
(of  course,  now  the  nondimensionalization  scheme  differs).  Using  the  transformations 
given  in  Sections  14  and  17,  we  replace  co  by  kco  and  d'r/d^  by  2a*d't,/dy  in  those 
relations.  Also,  +  cr$(2)  is  replaced  by  as  given  in  (199)  and  ,  f0  , 

f,  and  f2  are  given  by  (186),  (196),  and  (200),  respectively.  Under  this  procedure, 
tlie  relation  (102)  for  the  axial  component  of  the  velocity  perturbation  becomes,  taking 
note  of  the  form  of  (34),  as  follows: 
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Realizing  that  y  is  of  the  order  of  the  square  root  of  M  (the  Mach  number  along 
the  nozzle)  we  see  that  the  axial  component  of  the  perturbation  velocity  goes  to  zero 
as  (M)*  .  This  component  of  velocity  is  primarily  of  an  acoustic  nature  with  vorticity 
and  entropy  effects  of  the  order  of  crM2  and  CjM2  .  As  discussed  in  Reference  12, 
whenever  the  entropy  and  vorticity  waves  are  generated  by  the  combustion  process  they 
are  typically  of  0(M)  .  In  that  case,  cr  and  Cj  would  therefore  by  of  O(M)  , 
causing  the  effect  of  vorticity  and  entropy  waves  upon  the  axial  velocity  component 
to  be  of  0(M3)  . 
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The  relations  (102)  for  the  transverse  components  become 


(44)  shows  that  the  radial  velocity  perturbation  v'  behaves  as  pqrV  and  the 
tangential  velocity  perturbation  w'  behaves  as  (l/r)W  =  (l/r)V  .  From  the  steady- 
state  continuity  equation,  we  have  that  pqr*  is  a  constant  along  the  nozzle.  So, 
combining  rw/r  with  the  radial  dependence,  we  obtain  tha';  botlj  v'  and  wf  have 
axial  dependencies  which  behave  as  V/rw  or,  therefore,  r&  VM*  at  low  Mach  numbers. 
So,  each  of  the  two  transverse  velocity  components  are  primarily  affected  by  the 
entropy  and  vorticity  waves  rather  than  the  acoustic  waves.  The  rotational  effects 
have  both  an  undulatory  and  a  non-undulatory  portion.  These  effects  appear  to  be 
O(Mir)  while  the  acoustic  effect  is  0(M)  . 


The  thermodynamic  properties  now  remain  to  be  determined.  Following  the  same  pro¬ 
cedure  as  above,  (102)  gives  the  following  relation  for  the  axial  dependence  of  the 
pressure  perturbation: 
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The  pressure  perturbation  has  an  acoustic  behavior  which  is  of  0(Ml)  .  To  0(M3) 
there  is  no  undulatory  effect  due  to  the  entropy  and  vorticity  waves.  There  is, 
however,  a  non-undulatory  effect  of  0(a)  due  to  the  entropy  waves. 


Finally,  we  find  that  the  entropy  and  density  are  given  as 
S 

PR 


It  is  seen  that  the  entropy  wave  has  a  constant  amplitude  and  the  density  has  an 
undulating  effect  (due  to  entropy)  which  is  of  0(a)  . 

The  results  of  this  section  agree  qualitatively  with  the  results  of  Section  16. 
Figures  18,  19,  and  20  show  that  pressure  has  the  slightest  effect  due  to  entropy  and 
vorticity,  axial  velocity  has  a  small  effect,  and  transverse  velocities  have  a  large 
effect. 

It  should  be  noted  that  there  is  a  slight  inconsistency  throughout  the  last  two 
sections.  The  steady-state  density  was  considered  as  a  constant  of  0(1)  but, 
actually,  compressibility  causes  variations  of  0(y4)  .  This  means  that  any  terns 
in  the  previous  results  which  are  of  0(y“)  (or  0(M2))  higher  than  the  leading  terms 
cannot  be  considered  as  quantitatively  accurate.  The  reason,  however,  that  the 
analyses  were  continued  to  higher  orders,  was  to  establish  the  orders  of  magnitude 
of  the  acoustic,  entropy,  and  vorticity  effects.  *hile  the  steady-state  compres¬ 
sibility  would  modify  the  solution  to  higher  orders,  it  would  not  change  the  orders 
of  magnitudes  of  these  effects.  In  principle,  the  steady-state  density  variation 
could  have  been  considered  but  the  simplicity  of  (177)  would  have  been  destroyed. 

The  solutions  contain  several  groups  of  parameters:  (C?Beld)  ,  (X/B  e'^c+d^)  , 
(IoDkCj  +  a/2)  ,  (ak)  ,  (aX3)  ,  (a\j  +  CjXj)  ,  and  (\„/k)  .  The  first  two  parameters 
plus  a  and  Cj  are  given  by  initial  conditions  on  the  problem.  Then  the  other 
five  parameters  may  be  calculated.  Actually,  one  of  the  parameters  should  be  con¬ 
sidered  as  unitv  and  the  other  parameters  should  be  referenced  tc  it  since  the  problem 
is  linear. 
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Instead  of  calculating  the  asymptotes  for  and  as  described  in  Section 

17,  one  could  take  the  results  for  the  particular  solutions  of  (172)  and  substitute 
them  into  the  defining  relations  following  (162). 


19.  RESULTS  OF  THE  NOZZLE  ADMITTANCE  CALCULATIONS 
AND  THEIR  APPLICATIONS 

The  results  of  the  nozzle  admittance  calculations  are  presented  in  Tables  III  through 
LXXX.  The  values  of  the  real  and  imaginary  parts  of  the  coefficients  G  ,  8  ,  C  , 
a  ,  and  £  are  given  for  various  values  of  the  nozzle  entrance  Mach  number  M  ,  the 
angular  frequency  o)  ,  and  the  eigenvalue  s^j.  .  As  previously  noted,  the  frequency 
oi  must  be  multiplied  by  the  square  root  of  the  contraction  ratio  in  order  to  obtain 
(*)c  (the  frequency  nondimensionalized  by  the  chamber  radius).  The  coefficients  were 
calculated  at  constant  intervals  of  the  variable  0  which  do  not  correspond  to  con¬ 
stant  intervals  of  the  Mach  number. 

Of  course,  if  one  were  interested  in  values  of  s^  ,  co  ,  and/or  M  other  than 
those  presented  it  would  be  necessary  to  interpolate  or  extrapolate.  This  can  readily 
be  done,  will  be  shown  in  some  of  the  examples  given  later  in  this  section. 

The  most  important  admittance  coefficients  in  the  combustion  instability  application 
are  a  and  8  .  The  coefficient  a  depends  solely  upon  the  function  £  andtherefore 
does  not  contain  the  effects  of  vorticity  and  entropy  which  are  contained  in  , 

,  and  f3  .  8  is  a  combination  of  G  and  8  and  therefore  contains  neither 

nor  f3  .  Furthermore,  G  and  8  combine  in  such  a  way  so  as  to  suppress  the 
effect  of  £(i)  at  low  Mach  numbers.  So  8  contains  only  slight  effects  of  the 
vorticity  and  no  effects  of  the  entropy  waves.  For  these  reasons  a  and  8  do  not 
behave  as  widly  as  the  coefficients  G  ,  8  ,  and  C  and  interpolations,  are  more 
accurate  for  a  and  £  than  for  G  ,  8  ,  and  C  . 

Certain  relationships  have  been  developed  which  provide  necessary  and  sufficient 
conditions  for  the  neutral  stability  of  small  perturbations  in  a  rocket  combustion 
chamber*.  The  derivations  of  these  relationships  may  be  found  in  References  1,  3  and 
12.  A  typical  relationship  for  purely  transverse  oscillations  which  is  found  in 
Reference  12  is  as  followst: 
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(205) 


where  is  the  axial  dependence  of  the  energy  release  per  unit  volume,  L  is  the 

combustion  chamber  length,  ^  is  the  gas  velocity  to  stagnation  speed  of  sound  ratio 
at  the  nozzle  entrance,  and  cot  is  the  difference  between  the  angular  frequency  a>c 
and  the  eigenvalue  sl/h  .  A  similar  relationship  for  longitudinal  oscillations  which 
employs  a  instead  of  £  is  also  developed  in  Reference  12. 

*  Presently,  the  science  of  combustion  instability  has  practical  use  only  for  the  admittance 
coefficients  a  and  8  ;  G  and  8  have  been  used  only  in  the  combined  form  given  as  8 
and  C  have  been  neglected.  Here  a  more  general  viewpoint  was  adopted,  so  that  all  of  the 
admittance  coefficients  have  been  calculated  at  little  extra  expense. 

♦  Actually,  in  the  reference,  another  term  due  to  droplet  drag  has  been  Included.  However,  it 
is  not  important  for  the  purpose  of  the  example. 


79 


In  (205)  the  tern  on  the  left-hand  side  represents  the  integral  effect  over  the 
chamber  of  the  perturbation  of  energy  release  divided  by  the  pressure  perturbation. 

In  general,  it  will  be  a  complex  number  and  have  a  destabilizing  effect.  On  the  right- 
hand  side,  (y  +  1 )/y  has  a  stabilizing  effect  and  is  due  to  various  terms  in  the 
differential  equations  of  the  gas  motion.  The  second  term  on  the  right-hand  side 
represents  the  effect  of  the  nozzle  upon  the  stability  and  upon  the  frequency.  Whenever 
the  real  part  of  £  is  positive  it  has  a  destabilizing  effect,  while  whenever  the  real 
part  of  £  is  negative  it  has  a  stabilizing  effect.  The  last  term  represents  the 
change  in  frequency  from  the  acoustic  frequency.  It  has  no  effect  upon  the  stability 
since  it  is  imaginary. 

Now,  suppose  for  example  that  one  wishes  to  determine  the  stability  characteristics 
of  a  combustion  chamber  with  a  nozzle  of  the  shape  employed  in  the  calculations  (30° 
half-angle  in  the  conical  section  and  wall  radius  of  curvature  of  the  throat  equal  to 
the  throat  diameter).  The  nozzle  entrance  Mach  number  in  this  example  is  0.  113  and 
then  the  contraction  ratio  is  5.28.  The  first  tangential  mode  is  examined  jo  that 
s^jj  =  1.84  and,  taking  the  frequency  as  coc  =  1.72  ,  one  has  w  =  co,/(5.  28) 1  =  0.750  . 

In  order  to  find  the  values  of  £r  and  Si  ,  interpolations  of  the  results  presented 
in  the  tables  must  be  made.  The  following  summarizes  the  information  extracted  from 
the  tables  of  results  for  the  purpose  of  this  example: 


CO 

Si/h 

M 

Sr 

Si 

0.5 

1 

0.099 

0.074 

0.  114 

1.0 

1 

0.099 

-1.080 

-1.538 

0.5 

2 

0.099 

0.223 

1. 132 

1.0 

2 

0.099 

0.097 

0.063 

0.5 

1 

0.  124 

0. 100 

0.216 

1.0 

1 

0.  124 

-0.  383 

- 1. 008 

0.5 

2 

0.  124 

0.390 

1.319 

1.0 

2 

0.  124 

0.  130 

0.  230 

By  a  linear  interpolation  over  the  frequency  parameter  we  obtain,  for  oj  -  0.75  , 


Sj/Jj 

M 

£r 

£i 

1 

0.099 

-0.503 

-0.712 

2 

0.099 

0. 160 

0.598 

1 

0. 124 

-0. 142 

-0.  396 

2 

0.  124 

0.  260 

0.  775 

By  another  linear  interpolation  over  the  eigenvalue  parameter,  we  obtain  the 
following,  for  co  =  0.75  and  s^  =  1.84  , 


M 

Cl 

0.  099 

0. 054 

0.  388 

0. 124 

0.  195 

0.588 
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Now,  finally  interpolating  over  the  Mach  number,  we  have  for  co=  0.75  ,  =  1.84  , 

and  M  =  0.  113 

£r  =  0.133;  =  0.500. 

For  the  purpose  of  comparison,  the  equations  discussed  in  Section  14  were  integrated 
and  the  admittance  coefficients  were  calculated  for  the  values  w  =  0.75  and 
s^  =  1.84  .  With  M  =  0.113  that  exact  solution  was 

£r  =  0.136;  £t  =  0.481. 


This  compares  favorably  with  the  results  obtained  by  linear  interpolations;  there  is 
2%  error  in  £r  and  4%  error  in  £j  . 


There  are  many  cases  where  significant  errors  are  introduced  by  the  linear  inter¬ 
polation.  Often  the  situation  may  be  improved  by  employing  a  more  sophisticated  inter¬ 
polation  scheme.  For  example,  if  the  value  of  the  admittance  coefficient  at  three 
values  of  a  parameter  were  known,  a  parabola  could  be  fitted  to  determine  the  admittance 
coefficient  at  neighboring  values  of  the  parameters.  If  f  represents  the  particular 
admittance  coefficient  and  x  represents  a  particular  parameter,  the  following  formula 
results  from  such  a  parabolic  fit  to  the  points  xQ  ,  xx  ,  and  x2 


f(x) 


(x  -  xt)(x  -  x2) 
(xq  -  Xj)  (x0  -  x2) 


*<*.) 
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The  use  of  (206)  shall  be  demonstrated  in  an  example  where  interpolations  are  per¬ 
formed  only  over  the  svh  parameter.  Of  course,  it  may  be  used  successively  to  inter¬ 
polate  over  as  many  parameters  as  desired.  Consider  coc  =  2. 36  ,  s„h  =  1.84  ,  and 
M  =  0.250  .  The  contraction  ratio  is  2.45  so  th?t  wr  1. 50  .  The  following  summarizes 
the  pertinent  results  from  the  tables; 


sv>h 

Sr 

£i 

1 

-0.601 

-0.625 

2 

0. 116 

-0.060 

3 

0.  406 

0.689 

Using  the  results  for 
svh  =  1.84  . 


svh  =■ 


1 


and 


svb 


2  .  a  linear  interpolation  gives,  for 


0.001  ; 


-0. 154  . 
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On  the  other  hand.  (206)  may  be  used  with  x  representing  svb  and  f  representing 
£r  at  first  and  secondly  £j  .  Then  x0  =  1  .  x^  =  2  .  and  x2  =  3  .  The  result  is 
that,  for  s^  =  1.84  . 

£r  =  0.030;  8j  =  -0.163  . 

An  exact  solution  for  o>  =  1.50  .  s^h  =  1.84  ,  and  M  =  0.250  yields 
£r  =  0.036  ;  =  -0.198  . 

So  it  is  seen  that  the  more  sophisticated  interpolation  scheme  is  more  accurate  in 
this  case,  although  about  20%  error  -.till  exists*. 

In  both  cases  considered  above  £r  is  positive  so  that  the  nozzle  has  a  de¬ 
stabilizing  effect.  For  low  M  ,  Qi  M  ,  so  that  in  the  first  case  ^r/Y%  =  1-00 
while  in  the  second  case  it  equals  0. 120.  For  y  -  1. 2  ,  (y  +  *)/y  =  1.83  in  (205), 
so  that  the  destabilizing  effect  is  significant  by  comparison  in  the  first  case 
although  small  in  the  second  case. 

Suppose  now  that  we  did  not  have  a  30°  but  instead  a  15°  nozzle.  The  scaling 
factor  0  introduced  in' Section  12  is  chosen  as  (tan  15°)/(tan  30°)  =  0.465  .  In 
general  the  scaled  nozzle  will  not  be  identical  to  the  actual  nozzle  in  the  throat 
portion,  although  it  has  been  made  identical  in  the  conical  portion.  So.  while 
there  is  some  error  here  in  replacing  the  actual  nozzle  by  the  scaled  nozzle,  this 
error  is  expected  to  be  small  whenever  the  conical  portion  of  the  nozzle  is  con¬ 
siderably  larger  than  the  throat  portion. 

In  the  present  example  consider  the  throat  wall  radius  of  curvature  R  =  2  , 
coc  =  2. 97  .  s„h  =  1.0  ,  and  the  entrance  Mach  number  M  =  0. 152  .  Then  the  contraction 
ratio  is  3.94  and  cu  =  1.50  .  cj/0  equals  3.22  and  svb/0  =  2. 15  .  Since 

<3(«  .  w .  s„b)  =  ,  co/0,  s vXJ(5) 

and 

S(M  .  w ,  s„h)  =  /38ref  (M  ,  oo/0 ,  s„h//3)  , 

it  follows  that 

£(M  ,  co ,  s^jj)  =  £ref  (M  ,  oi/0  ,  s „h//3)  . 

It  remains  now  tc  determine  £ref  for  X  =  0. 152  ,  co  =  3.  22  ,  and  s^j,  =  2. 15  . 

From  the  tables  of  results,  one  obtains  the  following  information  for  M  =  0.  152  : 


Cl) 

Sr 

8i 

3.0 

2 

-0.  894 

-0.  134 

3.5 

2 

-0.982 

-0.  145 

3.0 

3 

-0.669 

-0.096 

3.5 

3 

-0.883 

-0.  193 

•  The  errors  are  considerably  saaller  when  the  more  proper  comparison  is  aade  with  the  differences 
between  the  values  of  6^.  (and  £1)  at  the  points  s^t  -  1.  2.  and  3  . 
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By  aeans  of  successive  linear  interpolations  over  the  frequency  and  eigenvalue  para¬ 
meters,  it  is  found  that,  for  M  =  0. 152  .  cj  =  3.  22  .  and  s„h  =  2. 15  . 

8r  =  -0.908;  =  -0.163. 


An  exact  calculation  for  these  values  of  cu  and  svh  with  8l  =  15°  and  R  =  2 
yields 

£r  =  -1.050;  Sj  =  -0.129, 


which  favorably  agrees  with  the  approximate  results  obtained  by  using  the  scaled 
nozzle.  The  difference  results  because  the  throat  wall  radius  of  curvature  R  for 
the  scaled  nozzle  is  greater  than  that  for  the  reference  nozzle  by  a  factor 
(l//3)z  =  4. 62  ,  while  the  actual  and  reference  nozzles  have  the  same  value  of  R  .  Of 
course,  if  the  actual  nozzle  had  an  identical  throat  contour  to  the  scaled  nozzle 
there  would  be  no  difference.  It  also  follows  that,  if  the  elliptical  arc  which  is 
tbe  generatrix  of  the  throat  wall  of  the  scaled  nozzle  (obtained  by  scaling  the  cir¬ 
cular  arc  of  the  reference  nozzle)  has  Identical  curvature  at  the  throat  to  the  curva¬ 
ture  of  the  actual  nozzle,  the  difference  in  the  results  would  be  quite  snail. 


In  certain  cases,  the  asymptotic  analysis  nay  provide  a  better  procedure  for  tbe 
deteralnation  of  the  adnittance  coefficients  than  Interpolation.  This  is  especially 
so  in  the  regions  of  wild  behavior  as  indicated  in  Figure  21.  There  the  irrotational 
adnittance  coefficient  a  is  plotted  versus  Mach  nuaber  along  the  nozzle.  It  is 
seen  that  the  asynptotic  solution  compares  favorably  with  the  exact  computer  solution. 
Wiereas  Interpolation  with  the  exact  results  over  Mach  nuxber  increments  of  0.025 
would  obviously  produce  serious  errors  in  tbe  region  of  wild  behavior,  the  asymptotic 
solution  wore  accurately  approximates  tbe  exact  solution. 


As  noted  in  Section  17,  the  constants  of  integration  in  the  asynptotic  solution 
must  be  determined  by  matching  the  exact  solution  to  the  aeyaptotic  solution  at  some 
point  in  the  low  Mach  nuaber  range.  In  Figure  21,  the  matching  point  is  M  =  0.099  . 
The  comparison  between  the  two  ablutions  is  most  favorable  when  the  matching  point  is 
in  or  near  the  region  of  wild  behavior.  Whenever  the  matching  point  is  away  from  this 
region  the  comparison  becomes  unfavorable. 

Furthermore,  tbe  comparison  for  l  is  better  than  for  and  .  Note 

£hat  a  depends  solely  upon  l  ,  while  the  other  admittance  coefficients  depend  upon 
■f(I)  and  sometimes  as  well.  Therefore  the  comparison  is  more  favorable  for 

a  thin  for  the  other  coefficients.  It  is  believed  that  the  unfavorable  comparisons 
for  £’(1)  and  £(2)  are  due  to  computational  difficulties  rather  than  theoretical 
difficulties.  At  this  time,  the  asymptotic  solutions  do  provide  a  substantial  amount 
of  insight  to  the  nature  of  tbe  oscillatory  flow  but  caution  must  be  used  whenever 
they  are  employed  to  actually  calculate  admittance  coefficients. 
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Flow 

Direction 


Nozzle  axis  — 


1  Geometry  of  convergent  portion  of  nozzle 


A(z,w,  $,h  )  =  Aref  -jf-  ) 

B(z,u?,$,h)  =  A  Br9f  (fiz,-^,^-  ) 
C(z  tWt$ )  -  Cref  (Pz,  -g-  t  p  ) 


Scaling  of  admittance  coefficients 
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Pi*. 7(c)  Real  part  of  pressure  adaittance  coefficient  versus  frequency 


Fig. 7(d)  Iaaginary  part  of  pressure  adaittance  coefficient  versus  frequency 


Frequency  cj 
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Fig.  10(a)  Real  part  of  irrotational  admittance  coefficient  versus  frequency 


Fig. 10(b)  Imaginary  part  of  irrotational  admittance  coefficient  versus  frequency 
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Fig. 11(a)  Real  part  of  combined  admittance  coefficient  versus  frequency 


Frequency  ou 


Fig.  11(b) 


Imaginary  part  of  combined  admittance  coefficient  versus  frequency 


Frequency  u 


Pig. 12(c)  Real  part  of  pressure  admittance  coefficient  versus  frequency:  Effect  of 

cone  angle 


Pig. 12(d)  Imaginary  part  of  pressure  admittance  coefficient  versus  frequency: 

Effect  of  cone  angle 
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Fig. 13(a)  Real  part  of  radial  velocity  admittance  coefficient  versus  frequency: 

Effect  of  throat  wall  curvature 
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Fig. 13(b)  Imaginary  part  of  radial  velocity  admittance  coefficient  versus  frequency: 

Effect  of  throat  wall  curvature 
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Pig. 13(c)  Real  part  of  radial  velocity  admittance  coefficient  versus  frequency: 

Effect  of  cone  angle 


Pig. 13(d)  Imaginary  part  of  radial  velocity  admittance  coefficient  versus  frequency: 

Effect  of  cone  angle 
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Pig.  14(a)  Real  part  cf  entropy  adaittance  coefficient  versus  frequency:  Effect  of 

throat  wall  curvature 


Fig. 14(b)  Imaginary  part  of  entropy  adaittance  coefficient  versus  frequency:  Effect 

of  throat  wall  curvature 
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Fig. 14(c)  Real  part  of  entropy  admittance  coefficient  versus  f-equency:  Effect  of 

cone  angle 


Frequency  u> 

Fig.  14(d)  Imaginary  part  of  entropy  admittance  coefficient  versus  frequency:  Effect 

of  cone  angle 
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Fig. 16(a)  Real  part  of  combined  admittance  coefficient  versus  fre 

throat  wall  curvature 


Fig. 16(b)  Imaginary  part  of  combined  admittance  coefficient  versui 

of  throat  wall  curvature 
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Fig. 16(c)  Real  part  of  combined  admittance  coefficient  versus  frequency:  Effect  of 

cone  angle 


Fig. 16(d)  Imaginary  part  of  combined  admittance  coefficient  versus  frequency:  Effect 

of  cone  angle 
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rfh  Nozzle  throat  Nozzle  axis 

— i - Nozzle  entrance 


Pig. 17  Nozzle  geometry  and  comparison  of  entrance  portions  of  approximate  and  actual 

nozzle  contours 


Axial  distance 


Fig.  18  Pressure  perturbation  versus  axial  distance  from  nozzle  throat 
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Fig.  19  Axial  velocity  perturbation  versus  axial  distance  from  nozzle  throat 
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Radial  velocity  perturbation  versus  axiel  distance  from  nozzle  throat 
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I  Irrotational  admittance  coefficient:  Comparison  between  exact  and  asymptotic 

solutions 


